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Abstract 



In this paper such Riemann metrics are estabhshed whose Laplace- 
Beltrami operators are identical to famihar Hamilton operators of ele- 
mentary particle systems. Such metrics are the natural positive definite 
invariant metrics defined on two-step nilpotent Lie groups. The corre- 
sponding wave and Schrodinger operators emerge in the Laplacians both 
of the static and solvable extensions of these nilpotent groups. The latter 
manifolds arc endowed with natural invariant indefinite metric of Lorentz 
signature. Thus, these new exact mathematical models provide a rela- 
tivistic theory for elementary particles. 

This theory establishes infinitely many non- equivalent models for which 
even classification is possible. The particle systems attached to them be- 
have exactly like their relatives introduced by the familiar standard model 
of elementary particle physics. Although there are strong connections be- 
tween the two theories, the new mathematical models were independently 
discovered. They appeared, originally, in those works of this author where 
isospectral Riemann manifolds with different local geometries have been 
constructed. Their strong connection to quantum theory, which seems 
to be completely unknown in the literature, was realized by some years 
later. 

One of the most important new features of this new theory is that the 
electromagnetic, the weak-nuclear, and the strong-nuclear forces emerge 
in a unified way. The main unifying idea is that these forces can be 
described by the eigenfunctions of the very same Laplacian such that 
the distinct forces emerge on distinct invariant subspaces of this common 
quantum operator. There are also other bonds realized which make the 
connection among these three fundamental forces much more strong. The 
missing fourth fundamental force, gravitation, is not discussed in this 
paper. 

On the solvable extensions the new models look like Fricdmann's ex- 
panding universe being adopted to the microscopic level. Like the macro- 
scopic one, also these microscopic expanding models obey Hubble's law. 
The microscopic models, however, offer much more complex structures 
with much more subtle explanations for some of those phenomena which 
have originally been clarified by the Friedmann model. For instance, it is 
rigorously established in this paper that, although the Riemann spaces in 
the mathematical models are not spatially isotropic in general, yet they 
are always spectrally isotropic. It follows that the radiation, which is 
experimentally known in cosmology, must also be isotropic. This state- 
ment actually contrasts the widely accepted view that the isotropic cos- 
mic radiation is the chief verification of the spatial-isotropy assumed by 
Friedmann in building up his theory. 
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1 Introduction 

By assuming the metric in the form ds"^ = a{t)'^ds^ — dt^, Friedmann (1922) 
built up his relativistic universe on two assumptions: the universe looks iden- 
tical in whichever direction we look, and this would also be true if we were 
observing it anywhere else. By the two differential equations imposed by 
the Einstein equation upon the function a(t), he showed then that, instead 
of being static, the universe is expanding. Without knowing Friedmann's 
theoretical prediction, this phenomenon was actually discovered by Edwin 
Hubble (1929) several years later. 

This phenomenon is pointed out in this paper in quite a different situation. 
Independently from Friedmann's theory, this one has been evolved from 
those works of the author where isospectral manifolds with different local 
geometries were constructed on two different types of manifolds, called Z- 
torus bundles (alias Z-crystals) and Z-ball bundles, respectively [Szl]-[Sz4]. 
These bundles are constructed by means of nilpotent Lie groups as well 
as their solvable extensions such that one considers tori resp. balls in the 
center (called also Z-space) of the nilpotent group. Both the nilpotent and 
extended groups are endowed with appropriate natural invariant metrics. 

Surprisingly enough, the Laplacians on the nilpotent groups (endowed al- 
ways with invariant positive definite Riemann metrics) are nothing but the 
familiar Hamilton operators corresponding to particle-antiparticle systems. 
On the two types of manifolds, the represented particles can be distinguished 
as follows. On Z-crystals, the Laplacian represents particles having no in- 
terior, where it actually appears as a Ginsburg-Landau-Zeeman operator 
of a system of electrons, positrons, and electron-positron-neutrinos. The 
particles represented by Z-ball bundles do have interior and the Laplacian 
decomposes into an exterior Ginsburg-Landau-Zeeman operator and an inte- 
rior spin operator by which the weak-force and the strong-force interactions 
can be described, respectively. These nuclear forces are very different from 
the electromagnetic force emerging in the Laplacian of Z-crystals. The weak 
nuclear force explains the beta decay, while the strong force keeps the parts 
of atomic nuclei together. Yet, the Z-crystal-Laplacian and the weak-force- 
Laplacian can be led back to the very same radial Ginsburg-Landau-Zeeman 
operator. This phenomenon is consistent with the Weinberg-Salam theory 
of beta decays, which unified the weak force with the electromagnetic force. 

There are two ways to introduce relativistic time on these models. The 
static model is constructed by the Cartesian product of the nilpotent group 
with M. The latter component becomes the time axis regarding the natural 
Lorenz-indefinite metric. According to the type of model being extended, the 
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Laplacian is the sum of Schrodinger and electron-positron- neutrino, resp., 
weak-nuclear and strong-nuclear wave operators. The last two wave op- 
erators appear for particles having interior. They are further decomposed 
into W- and Z-operators which are analogous to the electron-positron and 
electron-positron-neutrino wave operators. 

Relativistic time can be introduced also by solvable extension, which also 
increases the dimension of the nilpotent group by 1. The new axis, which is 
just a half-line M+, can also be used as time-axis for introducing a natural 
relativistic metric on these extensions. Contrary to the static case, in this 
way one defines expanding models obeying Hubble's law, furthermore, the 
Laplacian decomposes into expanding Schrodinger and electron-positron- 
neutrino, resp., weak and strong wave operators and the corresponding W- 
and Z-operators. 

It is a well known experimental fact that, even though the universe is 
expanding, there is no expansion measured on small scale level. Thus the 
question arises if the expanding solvable models describe real existing mi- 
croscopic world. Fortunately this question can positively be answered. In 
fact, despite of the expansion, the particles must not be expanding. The rea- 
son explaining this paradoxical phenomenon is that, defined by the angular 
momentum and spin operators, also these mathematical models correspond 
constant magnetic fields to the particles, and, due to the expansion, the 
change of these fields induces electromagnetic fields, which are completely 
radiated out from the system, keeping both the magnetic fields and the 
spectra of the particle-systems constant. There is also mathematically es- 
tablished in this paper that this radiation must be isotrpic, meaning that it 
is the same whichever direction is measured from. Thus the size and sev- 
eral other constants of particles must not be changing even according to the 
expanding solvable model. Actually this model gives a new explanation for 
the presence of an isotropic radiation in the space. 

Expanding on large scale but being stationary on small scale is a well 
known phenomenon, which, without the above explanation, could have been 
a major argument against the physical reality of the solvable extensions. 
According to physical experiments, although far distant clusters of galaxies 
move very rapidly away from us, the solar system is not expanding, nor is our 
galaxy or the cluster of galaxies to which it belongs. This stagnancy is even 
more apparent on microscopic level, where, for instance, the spectroscopic 
investigations of the light arriving from far distance galaxies confirm that 
the spectrum of hydrogen atom is the same today than it was macro billions 
of years ago. 
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The existence of isotropic background radiation is also well known which 
was measured, first, by Arno Penzias and Robert Wilson, in 1965. It is 
believed, today, that these radiations are travelling to us across most of 
the observable universe, thus the radiation isotropy proves that the universe 
must be the same in every direction, if only on a large scale. This phenome- 
non is considered as a remarkably accurate confirmation of Friedmann's first 
assumption. Our expanding model provides a much more subtle conclusion, 
however: This background radiation must be isotropic even if it arrives to 
us from very near distances. Moreover, it holds true also on non-isotropic 
spaces. 

This exact mathematical model is not derived from the standard model of 
elementary particles, which theory is based on a non-Abelian gauge theory 
where the basic objects are Yang-Mills connections defined on principal fibre 
bundles having structure group SU{3). Contrary to these gauge theories, 
in our case all physical quantities are defined by invariant Riemann metrics 
living on nilpotent resp. solvable groups. It will also be pointed out that no 
regular gauge-group exist on these models regarding of which these objects 
are gauge invariant. Yet, there is a bridge built up between the two theories, 
which explains why the particles introduced by the two distinct models 
exhibit the very same physical features. This bridge can be regarded as a 
correspondence principle associating certain Riemann metrics to the Yang- 
Mills models of elementary particles. 

The key point about the new exact mathematical model is that the center 
of the nilpotent group makes room to describe also the rich "inner life" of 
particles, which is known both experimentally and by theories explaining 
these experimental facts. This "inner life" is displayed by the de Broglie 
waves which appear in a new form in this new situation such that they 
are written up in terms of the Fourier transforms performed only on the 
center of the nilpotent group. By this reason, they are called Z-Fourier 
transforms, which are defined on the two types of models accordingly. On Z- 
crystals, where there is no "inner life" , it is nothing but the discrete Z-Fourier 
transform defined by the Z-lattice by which the Z-crystals are introduced. 
On Z-ball bundles, however, in order to obey the boundary conditions, more 
complicated so called twisted Z-Fourier transforms are introduced. 

The action of the very same Laplacian appears quite differently on these 
different function spaces. On Z-crystals, where there are no Z-boundary 
conditions involved, the strong nuclear forces do not appear either. In this 
case, the eigenfunctions arise as eigenfunctions of Ginsburg-Landau-Zeeman 
operators. By this reason, they are called electromagnetic eigenfunctions. 
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On Z-crystals, the theory corresponds to quantum electrodynamics (QED), 
while on Z-ball bundles it relates to quantum chromodynamics (QCD). 

In fact, on Z-ball bundles, due to the the Z-boundary conditions, the 
Laplacian appears in a much more complex form exhibiting both the weak 
and strong forces. More precisely, the weak force eigenfunctions satisfying 
a given boundary condition are defined by the eigenfunctions of the exte- 
rior Ginsburg-Landau-Zeeman operator introduced above for particles hav- 
ing interior. Although there are numerous differences between this exterior 
Ginsburg-Landau-Zeeman operator and the original GLZ-operator defined 
on Z-crystals, they are both reduced to the very same radial operator acting 
on radial functions. As a result, from the point of view of the elements of 
the spectrum, they are the same operators. This is the mathematical cer- 
tification of the Weinberg-Salam theory which unified the weak interaction 
with the electromagnetic force. The strong force eigenfunctions are defined 
by the eigenfunctions of the inner spin operator. All these forces reveal the 
very same strange properties which are described in QCD. This is how a 
unified theory for the three: 1.) electromagnetic, 2.) weak, and 3.) strong 
nuclear forces is established in this paper. The only elementary force missing 
from this list is the gravitation, which is not discussed in this paper. 

This very complex physical-mathematical theory can clearly be evolved 
just gradually. In order to understand the physical contents of the basic 
objects appearing in new forms in this new approach, first, those parts of 
the classical quantum theory are reviewed which are necessary to grasp 
these renewed versions of these basic concepts. Then, after introducing the 
basic mathematical objects on 2-step nilpotent Lie groups, several versions 
of the Z-Fourier transform will be studied. They are the basic tools both for 
introducing the de Broglie waves in a new explicit form and developing the 
theory unifying the three fundamental forces. Besides explicit eigenfunction 
computations, there is pointed out in this part that the Laplacian on Z- 
crystals is nothing but the Ginsburg-Landau-Zeeman operator of a system 
of electrons, positrons, and electron-positron-neutrinos. Furthermore, on 
the Z-ball models, it is the sum of the exterior Ginsburg-Landau-Zeeman 
operator and the interior spin operator by which the strong force interaction 
can be established. 

Then, relativistic time is introduced and both static and expanding mod- 
els are established. The Laplacian on these space-time manifolds appears 
as the sum of wave operators belonging to the particles the system consists 
of. The paper is concluded by pointing out the spectral isotropy in the 
most general situations. Since the Riemann metrics attached to the parti- 
cle systems are not isotropic in general, this statement points out a major 
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difference between our model and Friedmann's cosmological model where 
the isotropy of the space is one of his two assumptions. Our statement says 
that radiation isotropy holds true also on non-isotropic spaces and the two 
isotropy-concepts are by no means equivalent. 



2 Basics of classical quantum theory. 

In this sections three topics of classical quantum theory are reviewed. The 
first one describes the elements of de Broglie's theory associating waves to 
particles. The second resp. third ones are surveys on meson theory resp. 
Ginsburg-Landau-Zeeman and Schrodinger operators of charged particles. 

2.1 Wave-particle association. 

In quantum theory, a particle with energy E and momentum p is associated 
with a wave, ^e'^^'^~^*\ where K = (27r/A)n is the wave vector and n is the 
wave normal. These quantities yield the following relativistically invariant 
relations. 

For light quanta the most familiar relations are 
(1) E = huj , p = UK, 

where the length of the wave vector yields also the following equations: 



2 E ^_E^ 

2 ' "^^"^ If I ~ f — ' P ~ 2' 



(2) k=\K\=—, k = and \p\ = p = — , p 

c c c 



For a material particle of rest mass m, the fundamental relation is 

E 
c 



(3) — = \/p2 + m^c^ 



which can be established by the well known equations 



mc^ mv 



(4) E = , p - , 

of relativistic particle mechanics. 

The idea of de Broglie was that (1) should also be valid for a material 
particle such that (2) must be replaced by 



In this general setting, m = corresponds to the light. 
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In wave mechanics, de Broglie's most general wave packets are represented 
by the Fourier integral formula: 

(6) i^{Z,t) = j j I A{Ki,K2,K3)e'^^'''^^~'''^dK^dK2dK3, 

where uj is given by (5). In other words, a general wave appears as superpo- 
sition of the above plane waves. Instead of the familiar X, the vectors from 
the 3-space, M^, are denoted here by Z, indicating that the reformulated 
de Broglie waves will be introduced in the new theory in terms of the so 
called twisted Z-Fourier transform, which is performed just on the center 
(alias Z-space) of the nilpotent group. The X-space of a nilpotent group 
is complement to the Z-space and the integration in the formula of twisted 
Z-Fourier transform does not apply to the X-variable. It applies just to 
the Z-variable. The above denotation is intended to help to understand the 
renewed de Broglie waves more easily. 

The above wave function, ^p, satisfies the relativistic scalar wave equation: 

(7) {^'-^2§2)^iZ,t) = ^^l:{Z,t), 

which statement can be seen by substituting (6) into this equation. Then 
(5) implies (7), indeed. According to this equation, the wave function is 
an eigenfunction of the wave operator with eigenvalue m?c^/h?. By this 
observation we get that the spectrum of the wave operator is continuous 
and the multiplicity of each eigenvalue is infinity. 

The Fourier integral formula (6) converts differential operators to multi- 
plication operators. Namely we have: 

These correspondences together with (1) yield the translational key: 

(9) ' '^di^'' 

between the classical quantities p and E of classical mechanics and the 
operators of wave mechanics. 

In his lectures on physics [P] (Vol. 5, Wave mechanics, pages 3-4), 
Pauli describes the transition from the above relativistic theory to the non- 
relativistic approximation as follows. In mechanics, for v « c and p « 
mc, we have 

(10) - = Vp' + m^c^ ~ rnc{l + + ...) = -(mc^ + + . . . ). 
c 2 m^c^ c 2 m 
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From (5) we also obtain 

11 u = - = —- + —k' + ..., 

n n 2m 

where E = mc? + -Efcm and E^in = /2m. The non-relativistic wave 
(12) ^{Z,t) = j j I A{KuK2,K3)e'^^'''^^~'''^dK^dK2dK3, 
is defined in terms of 

which relates to the relativistic wave function by the formula: 

(14) i;{Z,t)=e-'-^'^[^iZ,t). 
Substitution into (7) yields then: 

, , „r, ~ m^c^ 7 im dib 1 d'^ib m^c^ ~ 

(15) vV + ^V' + 2-^-^^ = ^*, 

which is nothing but the non-relativistic wave equation: 
, , „9 ~ . 2m dih 1 d'^ib 



The imaginary coefficient ensures that there is no special direction in 
time. This equation is invariant under the transformations t — > — t and 
ijj ^ ijj* , whereby iptp*, where * means conjugation of complex numbers, 
remains unchanged (this conjugation will be denoted by ■0 later on). Ac- 
cording to quantum theory, the physically measurable quantity is not the 
wave functions ip or -i/;, but the probability densities resp. il^ip*. 



2.2 Meson theory. 



The above solutions of the wave equation strongly relate to the theory of 
nuclear forces and mesons. To review this field, we literally quote Hideki 
Yukawa's Nobel Lecture "Meson theory in its developments", delivered on 
December 12, 1949. Despite the fact that elementary particle physics went 
through enormous developments since 1949, this review has been chosen 
because it is highly suggestive regarding the physical interpretations of the 
mathematical models developed in this paper: 

"The meson theory started from the extension of the concept of the field 
of force so as to include the nuclear forces in addition to the gravitational 
and electromagnetic forces. The necessity of introduction of specific nuclear 
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forces, which could not be reduced to electromagnetic interactions between 
charged particles, was realized soon after the discovery of the neutron, which 
was to be bound strongly to the protons and other neutrons in the atomic 
nucleus. As pointed out by Wigner^, specific nuclear forces between two 
nucleons, each of which can be either in the neutron state or the proton 
state, must have a very short range of the order of 10-13 cm, in order to 
account for the rapid increase of the binding energy from the deuteron to the 
alphaparticle. The binding energies of nuclei heavier than the alpha-particle 
do not increase as rapidly as if they were proportional to the square of the 
mass number A, i.e. the number of nucleons in each nucleus, but they are 
in fact approximately proportional to A. This indicates that nuclear forces 
are saturated for some reason. Heisenberg^ suggested that this could be 
accounted for, if we assumed a force between a neutron and a proton, for 
instance, due to the exchange of the electron or, more generally, due to the 
exchange of the electric charge, as in the case of the chemical bond between a 
hydrogen atom and a proton. Soon afterwards, Fermi'^ developed a theory of 
beta-decay based on the hypothesis by Pauli, according to which a neutron, 
for instance, could decay into a proton, an electron, and a neutrino, which 
was supposed to be a very penetrating neutral particle with a very small 
mass. 

This gave rise, in turn, to the expectation that nuclear forces could be 
reduced to the exchange of a pair of an electron and a neutrino between 
two nucleons, just as electromagnetic forces were regarded as due to the 
exchange of photons between charged particles. It turned out, however, that 
the nuclear forces thus obtained was much too small^, because the betadecay 
was a very slow process compared with the supposed rapid exchange of the 
electric charge responsible for the actual nuclear forces. The idea of the 
meson field was introduced in 1935 in order to make up this gaps. Original 
assumptions of the meson theory were as follows: 

I. The nuclear forces are described by a scalar field U, which satisfies the 
wave equation 

in vacuum, where x is a constant with the dimension of reciprocal length. 
Thus, the static potential between two nucleons at a distance r is propor- 
tional to exp{—xr)/rh, the range of forces being given by 1/x. 

II. According to the general principle of quantum theory, the field U is 
inevitably accompanied by new particles or quanta, which have the mass 

(18) fi = — 
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and the spin 0, obeying Bose-Einstein statistics. The mass of these particles 
can be inferred from the range of nuclear forces. If we assume, for instance, 
X = 5x lO^^cm^, we obtain fi ~ 200me, where rue is the mass of the electron. 

III. In order to obtain exchange forces, we must assume that these mesons 
have the electric charge +e or — e, and that a positive (negative) meson is 
emitted (absorbed) when the nucleon jumps from the proton state to the 
neutron state, whereas a negative (positive) meson is emitted (absorbed) 
when the nucleon jumps from the neutron to the proton. Thus a neutron 
and a proton can interact with each other by exchanging mesons just as 
two charged particles interact by exchanging photons. In fact, we obtain an 
exchange force of Heisenberg type between the neutron and the proton of the 
correct magnitude, if we assume that the coupling constant g between the 
nucleon and the meson field, which has the same dimension as the elementary 
charge e, is a few times larger than e. 

However, the above simple theory was incomplete in various respects. 
For one thing, the exchange force thus obtained was repulsive for triplet 
S-state of the deuteron in contradiction to the experiment, and moreover we 
could not deduce the exchange force of Majorana type, which was necessary 
in order to account for the saturation of nuclear forces just at the alpha- 
particle. In order to remove these defects, more general types of meson 
fields including vector, pseudoscalar and pseudovector fields in addition to 
the scalar fields, were considered by various authors^. In particular, the 
vector field was investigated in detail, because it could give a combination 
of exchange forces of Heisenberg and Majorana types with correct signs and 
could further account for the anomalous magnetic moments of the neutron 
and the proton qualitatively. Furthermore, the vector theory predicted the 
existence of noncentral forces between a neutron and a proton, so that the 
deuteron might have the electric quadripole moment. However, the actual 
electric quadripole moment turned out to be positive in sign, whereas the 
vector theory anticipated the sign to be negative. The only meson field, 
which gives the correct signs both for nuclear forces and for the electric 
quadripole moment of the deuteron, was the pseudoscalar field''. There was, 
however, another feature of nuclear forces, which was to be accounted for as 
a consequence of the meson theory. Namely, the results of experiments on 
the scattering of protons by protons indicated that the type and magnitude 
of interaction between two protons were, at least approximately, the same as 
those between a neutron and a proton, apart from the Coulomb force. Now 
the interaction between two protons or two neutrons was obtained only if 
we took into account the terms proportional to g^, whereas that between a 
neutron and a proton was proportional to g'^, as long as we were considering 
charged mesons alone. Thus it seemed necessary to assume further: 
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IV. In addition to charged mesons, there are neutral mesons with the 
mass either exactly or approximately equal to that of charged mesons. They 
must also have the integer spin, obey Bose-Einstein statistics and interact 
with nucleons as strongly as charged mesons. This assumption obviously 
increased the number of arbitrary constants in meson theory, which could 
be so adjusted as to agree with a variety of experimental facts. These 
experimental facts could-not be restricted to those of nuclear physics in the 
narrow sense, but was to include those related to cosmic rays, because we 
expected that mesons could be created and annihilated due to the interaction 
of cosmic ray particles with energies much larger than cosmic rays in lOS?*^ 
was a great encouragement to further developments of meson theory. At that 
time, we came naturally to the conclusion that the mesons which constituted 
the main part of the hard component of cosmic rays at sea level was to be 
identified with the mesons which were responsible for nuclear force^. Indeed, 
cosmic ray mesons had the mass around 200me as predicted and moreover, 
there was the definite evidence for the spontaneous decay, which was the 
consequence of the following assumption of the original mesonn theory : 

V. Mesons interact also with light particles, i.e. electrons and neutrinos, 
just as they interact with nucleons, the only difference being the smallness 
of the coupling constant g in this case compared with g. Thus a positive 
(negative) meson can change spontaneously into a positive (negative) elec- 
tron and a neutrino, as pointed out first by Bhabha^O. The proper lifetime, 
i.e. the mean lifetime at rest, of the charged scalar meson, for example, is 
given by 



For the meson moving with velocity u, the lifetime increases by a factor 
— (i'/c)2 due to the well-known relativistic delay of the moving clock. 
Although the spontaneous decay and the velocity dependence of the lifetime 
of cosmic ray mesons were remarkably confirmed by various experiments^^, 
there was an undeniable discrepancy between theoretical and experimental 
values for the lifetime. The original intention of meson theory was to account 
for the beta-decay by combining the assumptions III and V together. How- 
ever, the coupling constant g, which was so adjusted as to give the correct 
result for the beta-decay, turned out to be too large in that it gave the life- 
time To of mesons of the order of 10~^sec, which was much smaller than the 
observed lifetime 2xW~^ sec. Moreover, there were indications, which were 
by no means in favour of the expectation that cosmic-ray mesons interacted 
strongly with nucleons. For example, the observed crosssection of scattering 
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of cosmic-ray mesons by nuclei was much smaller than that obtained theo- 
retically. Thus, already in 1941, the identification of the cosmic-ray meson 
with the meson, which was supposed to be responsible for nuclear forces, 
became doubtful. In fact, Tanikawa and Sakata^^ proposed in 1942 a new 
hypothesis as follows: The mesons which constitute the hard component of 
cosmic rays at sea level are not directly connected with nuclear forces, but 
are produced by the decay of heavier mesons which interacted strongly with 
nucleons. However, we had to wait for a few years before this two-meson 
hypothesis was confirmed, until 1947, when two very important facts were 
discovered. First, it was discovered by Italian physicistss that the negative 
mesons in cosmic rays, which were captured by lighter atoms, did not disap- 
pear instantly, but very often decayed into electrons in a mean time interval 
of the order of 10~^sec. This could be understood only if we supposed that 
ordinary mesons in cosmic rays interacted very weakly with nucleons. Soon 
afterwards, Powell and others^^ discovered two types of mesons in cosmic 
rays, the heavier mesons decaying in a very short time into lighter mesons. 
Just before the latter discovery, the two-meson hypothesis was proposed by 
Marshak and Bethe^^ independent of the Japanese physicists above men- 
tioned. In 1948 , mesons were created artificially in Berkeley and subse- 
quent experiments confirmed the general picture of two-meson theory. The 
fundamental assumptions are now^^ 

(i) The heavier mesons, i.e. n-mesons with the mass m^, about 280me 
interact strongly with nucleons and can decay into lighter mesons, i.e. vr- 
mesons and neutrinos with a lifetime of the order of 10~^sec; vr-mesons have 
integer spin (very probably spin 0) and obey Bose-Einstein statistics. They 
are responsible for, at least, a part of nuclear forces. In fact, the shape 
of nuclear potential at a distance of the order of h/mj^c or larger could be 
accounted for as due to the exchange of vr-mesons between nucleons. 

(ii) The lighter mesons, i.e. /i-mesons with the mass about 210me are the 
main constituent of the hard component of cosmic rays at sea level and can 
decay into electrons and neutrinos with the lifetime 2xl0~^ sec. They have 
very probably spin ^ and obey Fermi-Dirac statistics. As they interact only 
weakly with nucleons, they have nothing to do with nuclear forces. Now, if 
we accept the view that vr-mesons are the mesons that have been anticipated 
from the beginning, then we may expect the existence of neutral vr-mesons 
in addition to charged p-mesons. Such neutral mesons, which have integer 
spin and interact as strongly as charged mesons with nucleons, must be very 
unstable, because each of them can decay into two or three photons^ 8. In 
particular, a neutral meson with spin can decay into two photons and the 
lifetime is of the order of 10~^^sec or even less than that. Very recently. 
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it became clear that some of the experimental results obtained in Berke- 
ley could be accounted for consistently by considering that, in addition to 
charged n-mesons, neutral j^-mesons with the mass approximately equal to 
that of charged p-mesons were created by collisions of high-energy protons 
with atomic nuclei and that each of these neutral mesons decayed into two 
mesons with the lifetime of the order of W~^^sec or less^^. Thus, the neutral 
mesons must have spin 0. In this way, meson theory has changed a great deal 
during these fifteen years. Nevertheless, there remain still many questions 
unanswered. Among other things, we know very little about mesons heavier 
than TT- mesons. We do not know yet whether some of the heavier mesons 
are responsible for nuclear forces at very short distances. The present form 
of meson theory is not free from the divergence difficulties, although recent 
development of relativistic field theory has succeeded in removing some of 
them. We do not yet know whether the remaining divergence difficulties are 
due to our ignorance of the structure of elementary particles themselves^'^ . 
We shall probably have to go through another change of the theory, be- 
fore we shall be able to arrive at the complete understanding of the nuclear 
structure and of various phenomena, which will occur in high energy regions. 
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2.3 Zeeman and Schrodinger operators of electrons. 



In physics, the classical Zeeman operator of a charged particle is: 
(20) -l^^i^.v) - l^D, . +— -.(rr^ + y^) + eV. 



Originally, this operator is considered on the 3-space expressed in terms of 
3D Euclidean Laplacian and 3D magnetic dipole momentum operators. The 
latter operators are the first ones where a preliminary version of spin concept 
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appears in the history of physics. This is the so called exterior or orbiting 
spin associated with the 3D angular momentum 

(21) P = (Pl, P2, P3) = X P ^here 

(22) ^i = ^^^^3-^3P.) = }(Z2^-^3^).... 

Note that the above 2D operator keeps only component P3 of this angular 
momentum. The components of the complete 3D angular momentum obeys 
the commutation relations: 

(23) [Pi,P2]=iP3, [Pl,P3] = -iP2, [P2,P3] = iPl. 

In the mathematical models, the particles are orbiting in complex planes 
determined by the complex structures associated with the component an- 
gular momenta Pj, thus the 2D version plays more important role in this 
paper than the 3D version. The 2D operator is obtained from the 3D oper- 
ator by omitting Pi and P2 and by restricting the rest part onto the {x, y) 
plane. If Coulomb potential V is omitted, it is nothing but the Ginsburg- 
Landau-Zeeman operator of a charged particle orbiting on the (x, y)-plane 
in a constant magnetic field directed toward the z-axis. The magnetic dipole 
momentum operator, which is the term involving D^* ■= xdy — ydx and 
which is associated with the angular momentum operator hD^, commutes 
with the rest part, O, of the operator, therefore, splitting the spectral lines 
of O. The Zeeman effect is explained by this fine structure of the Zeeman 
operator. 

The Hamilton operator represents the total energy of a given physical 
system. More precisely, the eigenvalues of this operator are the discrete 
(quantized) energy values which can be assumed by the system. Thus, cor- 
respondence (9) implies Schrodinger's wave equation 

m -(^A,.„, - — . +^(.^ + - eV)^ = 

of an electron orbiting in the {x, y)-plane. 

As it well known, Schrodinger discovered first the relativistic equation 
which is a second order differential operator regarding the t-variable. By 
the time that Schrodinger came to publish this equation, it had already 
been independently rediscovered by O. Klein and W. Gordon. This is why 
it is usually called Klein-Gordon equation. Numerous problems had arise 
regarding this equation. Schrodinger became discouraged because it gave the 
wrong fine structure for hydrogen. Some month later he realized, however, 
that the non-relativistic approximation to his relativistic equation was of 
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value even if the relativistic equation was incorrect. This non-relativistic 
approximation is the famihar Schrodinger equation. 

Dirac also had great concerns about the Klein-Gordon equation. His 
main objection was that the probabilistic quantum theory based on this 
equation produced negative probabilities. Actually, the elimination of this 
problem led Dirac to the discovery of his relativistic electron equation. By 
this theory, however, proper probabilistic theory can be developed only on 
the relativistic space-time. This feature was strongly criticized by Pauli, 
according to whom such theory makes sense only on the space. 

The journey toward an understanding of the nature of spin and its rela- 
tionship to statistics has been taking place on one of the most difficult and 
exciting routes [B, Tom]. Although the Schrodinger wave equation gives ex- 
cellent agreement with experiment in predicting the frequencies of spectral 
lines, small discrepancies are found, which can be explained only by adding 
an intrinsic angular momentum to its usual orbital angular momentum of 
the electron that acts as if it came from a spinning solid body. The pioneers 
of developing this concept were Sommerfeld, Lande, and Pauli. They found 
that agreement with the Stern-Gerlach experiment proving the existence of 
Zeeman effect can be obtained by assuming that the magnitude of this addi- 
tional angular momentum was h/2. The magnetic moment needed to obtain 
agreement was, however, eh/2mc, which is exactly the same as that arising 
from an orbital angular moment of h. The gyromagnetic ratio, that is, the 
ratio of magnetic moment to angular momentum is therefore twice as great 
for electron spin as it is for orbital motion. 

Many efforts were made to connect this intrinsic angular momentum to 
an actual spin of the electron, considered it as a rigid body. In fact, the 
gyromagnetic ratio needed is exactly that which would be obtained if the 
electron consisted of a uniform spherical shell spinning about a definite axis. 
The systematic development of such a theory met, however, with such great 
difficulties that no one was able to carry it through to a definite conclusion. 
Somewhat later, Dirac derived his above mentioned relativistic wave equa- 
tion for the electron, in which the spin and charge were shown to be bound 
up in a way that can be understood only in connection with the requirements 
of relativistic invariance. In the non-relativistic limit, however, the electron 
still acts as if it had an intrinsic angular momentum of h/2. Prior to the 
Dirac equation, this non-relativistic theory of spin was originally developed 
by Pauli. 

Finally, we explain yet why the Coulomb operator does not appear in this 
paper and how can it be involved into the further investigations. The present 
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ignorance is mainly due to the fact that also the Hamilton operator (Lapla- 
cian) on nilpotent groups involves no Coulomb potential. There appear, 
instead, nuclear potentials like those Yukawa described in meson theory. An 
other major distinguishing feature is that this Laplacian (Hamilton oper- 
ator) includes also terms corresponding to the electron-positron-neutrino, 
which, by the standard model, is always something of a silent partner in 
an electron-positron-system, because, being electrically neutral, it ignores 
not only the nuclear force but also the electromagnetic force. Although the 
operator corresponding to this silent partner will be established by compu- 
tations evolved by Pauli to determine the non-relativistic approximation, all 
these operators appear in the new theory as relativistic operators complying 
with Einstein's equation of general relativity. 

The Coulomb force can be considered just later, after developing certain 
explicit spectral computations. This spectral theory includes also a spec- 
tral decomposition of the corresponding L^-function spaces such that the 
subspaces appearing in this decomposition are invariant with respect to the 
actions both of the Hamilton operator and the complex Heisenberg group 
representation. Also the latter representation is naturally inbuilt into these 
mathematical models. 

The complications about the Coulomb operator are due to the fact that 
these invariant subspaces (called also zones) are not invariant regarding the 
Coulomb's multiplicative operator [Sz5]- [Sz7]. In order to extend the the- 
ory also to electric fields, the Coulomb operator must be modified such that 
also this operator leaves the zones invariant. Such natural zonal Coulomb 
operator can be defined for a particular zone such that, for a given func- 
tion ip from the zone, function Vip is projected back to the zone. This 
modified Coulomb operator is the correct one which must be added to the 
Laplace operator on the nilpotent Lie group in order to have a relevant 
unified electro-magnetic particle theory. 

However, this modified Coulomb force is externally added and not nat- 
urally inbuilt into the Laplacian (Hamiltonian) of the nilpotent Lie group. 
In order to construct such Riemann manifolds whose Laplacian unifies the 
Ginsburg-Landau-Zeeman +neutrino-|- nuclear operators also with an ap- 
propriate Coulomb operator, the mathematical models must be further de- 
veloped such that, instead of nilpotent Lie groups, one considers general 
nilpotent-type Riemann manifolds and their solvable-type extensions. This 
generalization of the theory, which is similar to passing from special relativ- 
ity to the general one, will be the third step in developing this theory. 
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3 Launching the mathematical particle theory. 

The physical features of elementary particles are most conspicuously exhib- 
ited also by certain Riemann manifolds. A demonstration of this apparent 
physical content present in these abstract mathematical structures is, for 
instance, that the classical Hamilton and Schrodinger operators of elemen- 
tary particle systems appear as Laplace-Beltrami operators defined on these 
manifolds. Thus, these abstract structures are really deeply inbuilt into the 
very fabric of the physical world which can serve also as fundamental tools 
for building up a comprehensive unified quantum theory. 

Yet, this physical content of these particular mathematical structures has 
never been recognized in the literature so far. The ignorance is probably 
due to the fact that this new theory is not in a direct genetic relationship 
with fundamental theories such as the Gell-Mann~Ne'eman theory of quarks 
by which the standard model of elementary particles has been established. 
Neither is a direct genetic connection to the intensively studied super string 
theory, which, by many experts, is thought to be the first viable candidate 
ever for a unified quantum field theory of all of the elementary particles and 
their interactions which had been provisionally described by the standard 
model. Actually, the relationship between the two approaches to elementary 
particle physics is more precisely described by saying that there are both 
strong connections and substantial differences between the standard and our 
new models. In order to clearly explain the new features of the new model, 
we start with a brief review of the standard model. A review of string theory 
is, however, beyond the scope of this article. 

3.1 A rudimentary review of the standard model. 

The Gell-Mann~Ne'eman theory [GN] of quantum chromodynamics (QCD) 
grew out, in the 60's, from Yang-Mills' [YM] non-Abelian gauge theory 
where the gauge group was taken to be the SU (2) group of isotopic spin 
rotations, and the vector fields analogous to the photon field were inter- 
preted as the fields of strongly-interacting vector mesons of isotopic spin 
unity. QCD is also a non-Abelian gauge theory where, instead of SU{2), 
the symmetry group is SU{3). This model adequately described the rapidly 
growing number of elementary particles by grouping the known baryons and 
mesons in various irreducible representations of this gauge group. 

The most important difference between this and the new mathematical 
model is that the Yang-Mills theories do not introduce a (definite or indef- 
inite) Riemann metric with the help of which all other physical objects are 
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defined. They rather describe the interactions by Lagrange functions which 
contain more than a dozen arbitrary constants, including those yielding the 
various masses of the different kinds of particles. What makes this situation 
even worser is that all these important numbers are incalculable in principle. 
Under such circumstances the natural Hamilton and wave operators of par- 
ticles can not emerge as Laplacians on Riemann manifolds. This statement 
is the main attraction on the new models. 

To speak more mathematically, the fundamental structures for non-Abelian 
gauge theories are principal fibre bundles over Minkwski space with com- 
pact non-Abelian structure group SU{n) on which a potential is defined by 
a connection A, with components A^, in the Lie algebra su{n). The field is 
the curvature whose components are Fy^ = d^Ay — dyA^j, + [A^,Ay\. The 
most straightforward generalization of Maxwell's equations are the Yang- 
Mills equations dF = and d*F = 0, where d and d* are covariant deriva- 
tives. Gauge theories possess an infinite-dimensional symmetry group given 
by functions g : M ^ SU{n) and all physical or geometric properties are 
gauge invariant. 

To specify a physical theory the usual procedure is to define a Lagrangian. 
In quantum chromodynamics (QCD) such Lagrangian is to be chosen which 
is capable to portray the elementary particles in the following very rich 
complexity: The neutron and proton are composite made of quarks. There 
are quark of six type, or "flavors", the u, c, and t quarks having charge 2/3, 
and the d, s, and b quarks having charge —1/3 (these denotations are the first 
letters of words: up, down, charm, strange, top, and bottom). Quarks of 
each flavor come in three "colors" which furnish the defining representation 
3 of the SU (3) gauge group. 

Quarks have the remarkable property of being permanently trapped in- 
side "white" particles such as neutron, proton, baryon and meson. Baryon 
resp. meson are color-neutral bound states of three quarks resp. quarks 
and antiquarks. Neutron resp. proton are barions consisting an up- and 
two down- resp. one down- and two up-quarks. Thus the neutron has no 
charge while, in the same units in which the electron has an electric charge 
of —1, the proton has a charge of +1. The total charge of a white particle 
is always an integer number. Only the quarks confined inside of them can 
have non-integer charges. 

QCD can be regarded as the modern theory of strong nuclear forces hold- 
ing the quarks together. With no scalar fields, the most general renormaliz- 
able Lagrangian describing also these strong interactions can be put in the 
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form [W2] 

(25) C = --^F^^^-F^^, - Y,'^,,[x^^{d^ - \gA^^t^) + m„]V^„, 

n 

where ipnix) is a matter field, g is the strong couphng constant, ta are a 
complete set of generators of color SU (3) in the 3-representation (that is , 
Hermitian traceless 3x3 matrices with rows and columns labelled by the 
three quark colors), normalized so that Tr{tatp) = (5a^/2, and the subscript 
n labels quark flavors, with quark color indices suppressed. The first term 
is called matter Lagrangian density, the second one is the gauge field. 

Just as the electromagnetic force between electrons is generated by the 
virtual exchange of photons, so the quarks are bound to one another by a 
force that comes from the exchange of other quanta, called gluons because 
they glue the quarks together to make observable white objects. The gluons 
are flavor blind, paying no attention to flavor, however, they are very sen- 
sitive of color. They interact with color much as the photon interacts with 
electron charge. 

3.2 More specifics about the new abstract model. 

The above sketchily described objects are the most fundamental concepts 
in QCD. Their properties are established by the Lagrangians introduced 
there. In order to compare them, we review some more details about the 
new theory. The concepts introduced here will rigorously be establishment 
in the following sections. 

The mathematical structures on which the new theory is built on are 
2-step nilpotent Lie groups and their solvable extensions. Both type of 
manifolds are endowed with natural left invariant metrics. In this scheme, 
the nilpotent group plays the role of space, on which always positive defi- 
nite metric is considered. This choice is dictated also by the fact that the 
Hamilton operators of elementary particle systems emerge as the Laplacians 
of these invariant Riemann metrics. In order to ensure that these systems 
have positive energies, just positive definite metric can be chosen on these 
manifolds; for indefinite metrics the Laplacian never appears as the Hamil- 
ton operator of a particle system. 

Time can be introduced by adding new dimension to these nilpotent man- 
ifolds. This can be implemented either by a simple Cartesian product with 
the real line M, or by the solvable extensions of nilpotent groups. Both pro- 
cesses increase the dimension of the nilpotent groups by 1 and in both cases 
invariant indefinite metrics are defined such that the time-lines intersect the 
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nilpotent subgroup perpendicularly, furthermore, also {dt, dt) < holds. On 
these extended manifolds the Laplacian appears as the natural wave opera- 
tor (Schrodinger operator) attached to the particle systems. The difference 
between the two constructions is that the first one provides a static model, 
while the second one is an expanding model which yields the Hubble law of 
cosmology. 

Although both are relativistic, these space-time concepts are not quite 
the same than those developed in general relativity. Actually, Einstein's 
4D space-time concept has no room for exhibiting the rich "inner life" of 
particles which is attributed to them by meson theory, or, by the general 
standard model of elementary particle physics. This "inner life" can not be 
explained only by the properties of space-time. For instance, the symmetries 
underlying the electroweak theory are called internal symmetries, because 
one thinks of them as having to do with the intrinsic nature of the particles, 
rather than their position or motion. 

The abstract mathematical models, however, do make room for both the 
rich "inner life" and "exterior life" of particles. The main tool for exhibiting 
the inner physics is the center of the nilpotent Lie algebra, while the stage 
for the "exterior life" (that is, for the motion of particles) is the so called 
X-space denoted by X. This space is a complement of the center Z, which 
is called also Z-space. 

The space-like Z-space exhibits, actually, dualistic features. The primary 
meaning of vectors lying in the center is that they are the axes of angular mo- 
menta defined for the charged particles which are orbiting in complex planes 
in constant magnetic fields standing perpendicular to these complex planes. 
Actually, this axis-interpretation of vectors, Z, is developed in the following 
more subtle way: For any unit vector Z there is a complex structure, Jz 
acting on the X-space corresponded such that the particles are orbiting in 
the complex planes defined by Jz along the integral curves of vector fields 
defined by X — > Jz{X). As it is pointed out in the next section, this 2-step 
nilpotent Lie group is uniquely determined by the linear space, Jz, of skew 
endomorphisms Jz- Bijection Z ^ Jz provides a natural identification be- 
tween Z and Jz- More precisely, the group can be considered such that 
it is defined by a linear space of skew angular momentum endomorphisms 
acting on a Euclidean space, X, such that it is considered, primarily, as an 
abstract space Z which is identified with the endomorphism space, Jz, by 
the natural bijection Z — > Jz- Note that in this interpretation, the axis, Z, 
of the angular momentum, Jz, is separated from the complex plane where 
the actual orbiting is taking place. Anyhow, from this point of view, the 
Z- vectors exhibit space-like features. 
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Whereas, the constant magnetic field defined by the structure pins down 
a unique inertia system on which relations B = constant and E = holds. 
Thus a naturally defined individualistic inner time is given for each of these 
particles. This time can be synchronized, allowing to define also a com- 
mon time, T, which defines the time both on the static models and the 
solvable extensions. Prom this point of view, the center exhibits time-like 
features. This argument clarifies the contradiction between the angular- 
momentum-axis- and the customary time-axis-interpretation of the center 
of the Heisenberg groups. 

Although the concepts of relativity and quantum theory appear in new 
forms, they should be considered as refined versions of the original classical 
objects. Beside the above one, an other example for this claim is the new 
form by which de Broglie's waves are introduced on these groups. The most 
important new feature is that the Fourier transform is defined only on the 
center, Z = M.^, hy the following formula: 

(26) [ A{\X\,K)riz^'{Ku,X)zf{Ku,X)e'((^''^^-'^'UK, 

where, for a fixed complex basis B, complex coordinate system {zi{Ku, X)} 
on the X-space is defined regarding the complex structure Jk^ , for all unit 
Z- vector K^. This so called twisted Z-Fourier transform binds the Z-space 
and the X-space together by the polynomials Y\ -^f which depend both on 
the X- and K-variables. It appears also in several other alternative forms. 
Due to this complexity of the wave functions, the three main forces: the 
electromagnetic; the weak; and the strong forces of particle theory can be 
introduced in a unified way such that each of them can be expressed in terms 
of the weak forces. 

The main objects on these abstract structures are the Laplace operators 
considered both on nilpotent and solvable groups. They turn out to be the 
Hamilton resp. Schrodinger operators of the particle-systems represented 
by these metric groups. Due to the fact that these operators are Laplace 
operators on Riemann manifolds, the conservation of energy is automati- 
cally satisfied. In classical quantum theory the Hamilton function, which 
counts with the total energy of a system, is replaced by the Hamilton op- 
erator whose discrete eigenvalues are the quantized energy levels on which 
the system can exist. On the relativistic mathematical models the total 
energy is encoded into the Einstein tensor (stress-energy tensor) of the in- 
definite Riemannian metrics. In the quantum theory developed on these 
manifold, this stress-energy tensor is replaced by the Laplacians of these 
manifolds, which are, actually, the Hamilton resp. Schrodinger operators of 
the particle-systems represented by these models. In other words, this is a 
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correspondence principle associating the Laplacian resp. the eigenfunction- 
equations to the Einstein tensor resp. Einstein equation defined on these 
indefinite Riemann manifolds. 

Let it be emphasized again that this theory will be evolved gradually 
without adding any new objects to those defined mathematically on these 
abstract structures. The main focus is going to be to rediscover the most 
important physical features which are known by the standard model. Before 
starting this exploration, we describe, yet, a more definite bond between the 
two models. 



3.3 Correspondence principle bridging the two models. 

Correspondence principle associating 2-step nilpotent groups to SU{n)-mo- 
dels can also be introduced. It can be considered such that, to the Lagrange 
functions defined on the S'?7(n)-models, 2-step metric nilpotent Lie groups 
are corresponded. The combination of this correspondence principle with 
the above one associates the Laplacian of the metric group to the Lagrange 
function defined on the Yang-Mills model. This association explains why 
are the conclusions about the nature of electromagnetic, strong, and weak 
forces so similar on the two models. This bridge can be built up as follows. 

As it is explained above, the invariant Riemann metrics defined on 2- 
step nilpotent groups, modelling the particle systems in the new theory, 
can be defined for any linear space, Jz, of skew endomorphisms acting on 
the X-space. Thus, for a faithful representation, p, of su{n) C so{2n) in 
the Lie algebra of real orthogonal transformations acting on a Euclidean 
space X, one can define a natural 2-step nilpotent metric group by the 
endomorphism space Jz = p{su[n)) C so{X) whose X-space is X and the 
Z-space is the abstract linear space Z = p{su[n)). Actually, this is the 
maximal 2-step nilpotent metric group which can be corresponded to a Yang- 
Mills principal fibre bundle having structure group SU{n), and whose fibres 
consists of orthonormal frames of X on which the action of p{SU (n)) is one- 
fold transitive. Note that this group is still independent from the Yang-Mills 
connection A^ by which the Yang-Mills field (curvature), F^^, is defined. 

By the holonomy group Holp{A) C p{SU{n))p, defined at a fixed point p 
of the Minkowski space, groups depending on Yang-Mills fields can also be 
introduced. Even gauge-depending groups can be constructed, by sections 
cr : M — > M, where M denotes the total space of the fibre bundle and 
is lying in the fibre over the point p. In this case the endomorphism space 
is spanned by the skew endomorphisms A^{a{p)) considered for all p G M 
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and indices /x. This construction depends on sections cr(p). Since the gauge 
group is transitive on the set of these sections, this correspondence is not 
gauge invariant. 

These correspondences associate 2-step nilpotent metric groups also to the 
representations of the particular Lie algebras su{2) C so(4) resp. su{3) C 
so{6) by which the Yang-Mills- resp. Gell-Mann~Ne'eman-models are in- 
troduced. This association does not mean, however, the equivalence of the 
two theories. The X-space (exterior world) of the associated group is the 
Euclidean space X where the skew endomorphisms from p[su{n)) are act- 
ing. The Z-space is the abstract space Z = p[g) where g can be any of 
the subspaces of su{n) which were introduced above. In other words, for 
a Yang-Mills model, X-space X is already given and angular momentum 
endomorphism space is picked up in p[su{n)) in order to define the Z-spaces 
of the corresponded nilpotent groups. The space-time, including both the 
exterior and interior worlds, are defined by these spaces. Note that this con- 
struction completely ignores the Minkowski space which is the base-space for 
the Yang-Mills principle fibre bundle. These arguments also show, that the 
nilpotent groups corresponded to a fixed Yang-Mills model are not uniquely 
determined, they depend on the Z-space chosen on the YM-model. 

Whereas, on a YM-model, the exterior world is the Minkowski space 
over which the principal fibre bundle is defined and the interior world is 
defined there by p{su{n)). An other fundamental difference is that YM- 
models are based on submersion-theory, however, not this is the case with 
the new models. For instance, the Hamilton operators of particles never 
appear as sub-Laplacians defined on the X-space. Quite to the contrary, 
these Hamilton operators are acting on the total {X, Z)-space, binding the 
exterior and interior worlds together into an unbroken unity not character- 
istic for submersions. These are the most important roots explaining both 
the differences and similarities between the two models. These arguments 
also show that the new theory deals with much more general situations than 
those considered in Yang-Mills' resp. Gell-Mann~Ne'eman's SU{2)- resp. 
5'C/(3)-theories. It goes far beyond the SU{n) x • • • x 5?7(n)-theories. Since 
the Laplacians are natural Hamilton operators of elementary particles also 
in these most general situations, there is no reason to deny their involvement 
into particle theory. 

Differences arise also regarding the Maxwell theory of electromagnetism. 
The S'C/(n)-theories are non-Abelian gauge theories where the field is the 
curvature, F, of a Yang-Mills connection (potential) and the properties are 
gauge invariant regarding an infinite dimensional gauge group whose Lie 
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algebra consists of sn(n)-valued 1-forms, lo, satisfying div = F. The cus- 
tomary reference to this phenomenon is that only F is the physical object 
and 1-forms w do not have any physical significance, they are the results 
of mere mathematical constructions. This interpretation strongly contrasts 
the Aharanov-Bohm theory where these vector potentials do have physical 
meanings by which the effect bearing their names can be established [AB], 
[T]. Whereas, on the nilpotent groups and their solvable extensions the 
fundamental fields are the natural invariant Riemann metrics, g, by which 
all the other objects are defined. The classical Hamilton resp. Schrodinger 
operators emerge as Laplacians of metrics g. Curvature F corresponds to 
the field g{Jz{X),Y) in this interpretation. Since the basic objects are not 
invariant regarding their actions, the gauge-symmetries are not involved to 
these investigations. One-form uJziY) = {l/2)g{Jz{X), Y), defined over the 
points X £ JY, is the only element of the Lie algebra of the gauge-symmetry 
group which is admitted to these considerations. It defines the constant 
magnetic field and vanishing electric field which are associated with the 
orbiting spin and inner time T. 

This interpretation shows that this model breaks off, at some point, from 
the Maxwell theory whose greatest achievement was that it unified the, until 
1865, separately handled partial theories of magnetism and electricity. The 
Yang-Mills gauge theory is a generalization of this unified theory to vector 
valued fields and potentials. The nilpotent Lie group model approaches to 
the electromagnetic phenomena from a different angle. Since the potentials 
do have significance there, it stands closer to the Aharanov-Bohm theory 
than to the Maxwell- Yang-Mills gauge theory. Actually the full recovery 
of electromagnetism under this new circumstances will be not provided in 
this paper. Note that only the magnetic field has been appeared so far 
which is not associated with a non-trivial electric field. In other words, the 
magnetism and electricity emerge as being separated in this paper. The 
reunion of this temporarily separated couple can be established just at a 
later point after further developing this new theory. 

This bridge explains the great deal of properties which manifest similarly 
on these two models, however, turning from one model to the other one is 
not simple and the complete exploration of overlapping phenomena requires 
further extended investigations. In this paper the first order task is to firmly 
establish the point about the new model. Therefore, when several properties 
are named by the same name used in QCD, their definition strictly refer to 
our setting. Whereas, by these deliberately chosen names, the similarities 
between the two theories are indicated. 
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4 Two-step nilpotent Lie groups. 
4.1 Definitions and interpretations. 

A 2-step nilpotent metric Lie algebra, {M, (,)}, is defined on a real vector 
space endowed with a positive definite inner product. The name indicates 
that the center, Z, can be reached by a single application of the Lie bracket, 
thus its second application always results zero. The orthogonal complement 
of the center is denoted by X. Then the Lie bracket operates among these 
subspaces according to the following formulas: 

(27) [J\f,J\f]=Z , [J\f,Z]=0 , J\f = X®Z = R''xRK 

Spaces Z and X are called also Z- and X-space, respectively. 

Upto isometric isomorphisms, such a Lie algebra is uniquely determined 
by the linear space, Jz, of skew endomorphisms Jz '■ X ^ X defined for 
any Z £ Z hy the formula 

(28) {[X,Y],Z) = {Jz{X),Y)yZ eZ. 

This statement means that for an orthogonal direct sum, J\f = X (B Z = 
M.^ X R', of Euclidean spaces a non-degenerated linear map, J : ^ — > 
SE(X) , Z — > Jz, from the Z-space into the space of skew endomorphisms 
acting on the X-space, defines a 2-step nilpotent metric Lie algebra on M by 
(28). Furthermore, an other non-degenerated linear map J having the same 
range Jz = as J defines isometrically isomorphic Lie algebra. 

By means of the exponential map, also the group can be considered 
such that it is defined on M. That is, a point is denoted by {X, Z) on 
the group as well. Then, the group multiplication is given by the formula 
{X, Z){X', Z') = {X + X',Z + Z' + l[X.X']). Metric tensor, g, is defined 
by the left invariant extension of (, ) onto the group N. 

Endomorphisms Jz(-) will be associated with angular momenta. It must 
be pointed out, however, a major conceptual difference between the classical 
3D angular momentum, introduced in (22), and this new sort of angular 
momentum. For a fixed axis Z £ the endomorphism associated with 
the classical 3D angular momentum is defined with the help of the cross 
product X by the formula Jz '■ X ^ Z x X . That is, axis Z is lying in 
the same space, M^, where the endomorphism itself is acting. Linear map 
J : Z — > SE{X) on a 2-step nilpotent Lie group, however, separates axis 
Z € Z from the X-space where the endomorphism Jzi-) is acting. In other 
words, the latter endomorphism defines just orbiting of a position vector X 
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in the plane spanned by X and Jz{X), but the axis of orbiting is not in the 
X- but in the Z-space. 

In this respect, the Z-space is the abstract space of the axes associated 
with the angular momentum endomorphisms. According to this interpre- 
tation, for a fixed axis Z in the Z-space, a particle occupies a complex 
plane in the complex space defined by the complex structure Jz on the X- 
space. Abstract axis, Z, is considered as an "inner dial" associated with 
the particles, which is represented separately in the Z-space. This Z-space 
contributes new dimensions to the X-space which is considered as the inner- 
world supplemented to the exterior-world in order to have a natural stage 
for describing the inner physics of elementary particles. 

The above definition of 2-step nilpotent Lie groups by their endomorphism 
spaces Jz shows the large variety of these groups. For instance, if Z is an 
/-dimensional Lie algebra of a compact group and J : — > so{k) (which 
corresponds Jz G so{k) to Z £ Z) is its representation in a real orthogonal 
Lie algebra (that is, in the Lie algebra of skew-symmetric matrices) defined 
for X = M.^, then the system {Af = X (B Z, Jz} defined by orthogonal 
direct sum determines a unique 2-step nilpotent metric Lie algebra where 
the inner product on Z is defined by {Z,V) = —Tr{Jz o Jy). Thus, any 
faithful representation JJ : Z — > so{k) determines a unique two-step nilpotent 
metric Lie algebra. Since Lie algebras su{k/2) C so{k) used in non-Abelian 
gauge theories are of compact type, therefore, to any of their representations 
in orthogonal Lie algebras, one can associate a natural two-step nilpotent 
metric Lie algebra. This association is the natural bridge between a non- 
Abelian S'C/(n)-theory and the new theory developed in this paper. 

The 2-step nilpotent Lie groups form even a much larger class than those 
constructed above by orthogonal Lie algebra representations. In fact, in the 
most general situation, endomorphism space Jz is just a linear space defined 
by the range of a non-degenerated linear map I : Z ^ SE{X) which may 
not bear any kind of Lie algebra structure. However, such general groups 
can be embedded into those constructed by the compact Lie algebras, Jj. 
The smallest such Lie algebra is generated by the endomorphism space Jz 
by the Lie brackets. 

Very important particular 2-step nilpotent Lie groups are the Heisenberg- 
type groups, introduced by Kaplan [K], which are defined by endomorphism 
spaces Jz satisfying the Clifford condition j| = —z'^id, where z = \Z\ de- 
notes the length of the corresponding vector. These groups are attached 
to Clifford modules (representations of Clifford algebras). The well known 
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classification of these modules provides classification also for the Heisenberg- 
type groups. According to this classification, the X-space and the endomor- 
phisms appear in the following form: 

(29) X = {W'^^'^r X i^^^y ,Jz = {jzx...x jz) X i-jz X • • • X -jz), 

where / = dim(Z) and the endomorphisms jz act on the corresponding com- 
ponent, M''^'\ of this Cartesian product. The groups and the corresponding 
natural metrics are denoted by H^""'^^ and g^""'^^ respectively. Particularly 

important examples are the H-type groups H^'^''^\ where the 3-dimensional 
Z-space, M^, is considered as the space of imaginary quatrnions, further- 
more, action of jz on the space M^(^) = EI = of quaternoinic numbers is 
defined by left multiplications with Z. 

A brief account on the classification of Heisenberg type groups is as fol- 
lows. If / = dim{Jz) / 3 mod 4, then, upto equivalence, there exist exactly 
one irreducible H-type endomorphism space acting on a Euclidean space 
]R"i, where the dimensions n;, which depend just on Z, are described below. 
These endomorphism spaces are denoted by j/^"*- If i = 3 mod 4, then, 
upto equivalence, there exist exactly two non-equivalent irreducible H-type 
endomorphism spaces acting on M"'. They are denoted by and jf^'^^ 

respectively. They relate to each other by the relation J^^'^^ — 

The values n; corresponding to / = 8p, 8p + 1, . . . , 8p + 7 are 



(30) ni = 2^^ 2^^+"^ 2^P+^ 2^^+^ 2^^+^ 2"'^'+^ 2^^+^ 2'^^+'^ 

The reducible Clifford endomorphism spaces can be built up by these 
irreducible ones. They are denoted by Ji°^ resp. J^'^'^\ The corresponding 
Lie algebras are denoted by Tii""^ and Tti""'^^ respectively, which define the 

groups Hr^^ resp. H^"''^\ In the latter case, the X-space is defined by the 
(a + 6)-times product M'"' x • • • x M"' such that, on the last h component, the 

action of a Jz is defined by — —J^^''^\ and, on the first a components, 

the action is defined by J^z'^'^ ■ ^^^^ process should be applied 

only on the corresponding a-times product. 

One of the fundamental statements in this theory is that, in case of Z = 3 
mod 4, two groups H^""'^^ and Hi"" ''' ■* are isometrically isomorphic if and 
only if (a, b) = (a', b') upto an order. By a general statement, two metric 2- 
step nilpotent Lie groups with Lie algebras M = X QZ and M' = X' ®Z' are 
isometrically isomorphic if and only if there exist orthogonal transformations 
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A : X ^ X' and B : Z ^ Z' such that J^i^z) = A o o holds, for 

ah Z ^ Z. The isomorphic isometry between H^'^'^^ and ffj-^'""^ is defined 
hy A = id and B = —id. If (a, b) 7^ (a', b') (upto an order) then the 
corresponding groups are not isometrically isomorphic. 

In order to unify the two cases, denotations jj^'^^ = J^^^ and J^'^'^^ = 

—Jj^^ are used also in cases I 7^ 3 mod 4. One should keep in mind, however, 
that these endomorphism spaces are equivalent, implying that two groups 
H^'^'^^ and H^'^ ^ defined by them are isometrically isomorphic if and only 
[{ a + b = a' + b' holds. 

H-type groups can be characterized as being such particular 2-step metric 
nilpotent Lie groups on which the skew endomorphisms, Jz, for any fixed 
Z £ Z, have the same eigenvalues itzi. By polarization we have: 

(31) + JzM = -{Zi,Z2)Id, 

which is called Dirac's anticommutation equation. It implies that two en- 
domrphisms, Jzi and Jz2, with perpendicular axes, Zi _L Z2, anticommute 
with each other. Endomorphism spaces satisfying this weaker property de- 
fine more general, so called totally anticommutative 2-step nilpotent Lie 
groups on which the endomorphisms can have also properly distinct eigen- 
values. The classification of these more general groups is unknown in the 
literature. 

Let it be mentioned, yet, that groups constructed above by SM(2)-repre- 

sentations are exactly the groups h'^'^\ while those constructed by su(3)- 
representations are not even totally commutative spaces. Thus, they are 
not H-type groups either. It is also noteworthy, that the Cliffordian endo- 
morphism spaces J^^"'*^ do not form a Lie algebra in general. In fact, only 
the endomorphism spaces defined for I = 1,3,7 can form a Lie algebra. In 
the first two cases, pair (a, b) can be arbitrary, while in case / = 7 only 
cases = 1,6 = 0, or, a = 0,6 = 1 yield Lie algebras. It is an interesting 
question that which orthogonal Lie algebras can be generated by Cliffordian 
endomorphism spaces J/(a, 6). Let it be mentioned, yet, that Lie algebra 
sn(3) does not belong even to this category. 

4.2 Laplacian and curvature. 

Although most of the results of this paper extend to the most general 2- 
step metric nilpotent Lie groups, in what follows only H-type groups will be 
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considered. On these groups, the Laplacians appear in the following form: 

1 

(32) A = Ax + (1 + 4x2)A^ + 

a=l 

where D^* denotes directional derivatives along the vector fields X — > 
Ja{X) = Je^{X), furthermore, x = \X\ denotes the length of X-vectors. 

This formula can be established by the following explicit formulas. Con- 
sider orthonormal bases {Ei; . . . ; E^} and {ei; . . . ; e^} on the X- and Z- 
space respectively. The coordinate systems defined by them are denoted by 
{x^; . . . ; x'^} and {z^; . . . ; 2;'} respectively. Vectors extend into the 

left-invariant vector fields 

I I 

(33) Xi = di + - ^([X, Ei],ea)da = d^ + -Y,{J^ [X) , Ei)d^ 

Q=l Q=l 

and Zq = respectively, where di = d/dx^^ da = d/dz°' and J a = Je^- 

The covariant derivative acts on these invariant vector fields according to 
the following formulas. 

(34) VxX* = '^[X,X*] , VxZ = VzX = -]^Jz{X) , V^Z* = 0. 

The Laplacian, A, acting on functions can explicitly be established by 
substituting (33) and (34) into the following well known formula 

k I 

(35) A = ^ (Xf - Vx,X,) + J] (Z2 - Vz^Z,). 

i=l a=l 

These formulas allow to compute also the Riemannian curvature, R, on 
N explicitly. Then we find: 

(36) R{X,Y)X* = lj[x,Y]iXn - \jlY,x*]{X) + \j[x,x*]{Y); 
{37)R{X,Y)Z = -^[X, Jz{Y)] + j[Y, Jz{X)]; R{Zi, ^2)^3 = 0; 

(38) R{X,Z)Y = -^[X,Jz{Y)]; R{X, Z)Z* = -^JzJz^iX); 

(39) R{Z,Z*)X = --^Jz.JziX) + -^JzJz'iX), 

where X]X*]Y G X and Z; Z*; Zi; Z2; Z3 G Z are considered as the ele- 
ments of the Lie algebra N . The components of this tensor field on coordi- 
nate systems {x^; . . . ; x*^; z^; . . . ; z'} can be computed by formulas (33). 



32 



EXACT MATH-MODELS OF MICRO UNIVERSES 



5 Particles without interior. 

These Riemann manifolds were used, originally [Szl]-[Sz4], for isospectral- 
ity constructions in two completely different situations. In the first one, 
the Z-space is factorized by a Z-lattice, Tz, defined on the Z-space, which 
process results a torus bundle over the X-space. In the second case, Z-ball 
resp. Z-sphere bundles are considered by picking Z-balls resp. Z-spheres 
in the Z-space over the points of the X-space. It turns out that, apart 
from a constant term, the Laplacian on the Z-torus bundles, called also Z- 
crystals, is the same as the classical Ginsburg-Landau-Zeeman operator of 
an electron-positron system whose orbital angular momentum is expressed 
in terms of the endomorphisms Jz^ defined by the lattice points Zq, S Tz- 
More precisely, the Z-lattice defines a natural decomposition Wa of the 
L^-function space such that the components Wa are invariant under the ac- 
tion of the Laplacian, which, after restricting it onto a fixed Wa, appears as 
the Ginsburg-Landau-Zeeman operator whose orbital angular momentum is 
associated with the fixed endomorphism Jz^- 

It turns out, in the next chapters, that the constant term corresponds 
to neutrinos (massless particles with no charge) accompanying an electron- 
positron system. Thus, altogether, the Laplacian appears as the Hamilton 
operator of a system of electrons positrons and electron-positron-neutrinos 
which is formally the sum of a Ginsburg-Landau-Zeeman operator and a 
constant term. Names Ginsburg-Landau indicate that no Coulomb potential 
or any kind of electric forces are involved into this operator. Thus the 
forces which manifest themself in the eigenfunctions of the Laplacian, are 
not the complete electromagnetic forces, yet. However, when also these 
forces will be introduced, the eigenfunctions remain the same, the electric 
force contributes only to the magnitude of the eigenvalues. By this reason, 
the forces associated with these models are called electromagnetic forces. 
Since the Z-lattices consist of points and intrinsic physics is exhibited on 
the Z-space, particles represented by Z-crystals are considered as point-like 
particles having no insides. The theory developed for them is in the strongest 
connection with quantum electrodynamics (QED). 

In the second case both the Laplacian and the angular momentum oper- 
ator appear in much more complex forms. Contrary to the first case, the 
angular momentum operator is not associated with a fixed Z-vector, but 
it represents spinning about each Z-vector. Beside the orbiting spin, also 
natural inner operators emerge which can be associated both with weak and 
strong nuclear forces. The particles represented by these models do have 
inside on which the intrinsic physics described in QCD is exhibited on a full 
scale. An openly admitted purpose in the next sections is to give a unified 
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theory for the 3 forces: 1.) electromagnetic- 2.) weak-nuclear- 3.) strong- 
nuclear- forces. However, some concepts such as Dirac's spin operator will 
be introduced in a subsequent paper into this theory. 

5.1 Z-crystals modelling Ginsburg-Landau-Zeeman operators. 

The Z-torus bundles are defined by a factorization, T\H, of the nilpotent 
group -ff by a Z-lattice, T = {Z^}, which is defined only on Z and not on 
the whole {X, Z)-space. Such a factorization defines a Z-torus bundle over 
the X-space. The natural Z-Fourier decomposition, := W-y, of the 
function space belonging to this bundle is defined such that subspace is 
spanned by functions of the form 

(40) ^'^(X, Z) = ^{X)e'^'''^^'"^l 

Each is invariant under the action of A, more precisely we have: 

(41) A^^(X,Z) = (<^V)(^)e^^'^^^'^\ where 

(42) <^ = Ax + 2^11)^ .-4^222(1 + ^x2). 

In terms of parameter = ttz^, this operator is written in the form <^ = 
Ax + 2iDf^ • — /i^x^ — 4//2. Although it is defined in terms of the X-variable, 
this operator is not a sub-Laplacian resulted by a submersion. It rather 
is the restriction of the total Laplacian onto the invariant subspace W-y. 
Actually, the Z-space is represented by the constant /i and operator 
A characteristic feature of this restricted operator is that it involves only a 
single endomorphism, Jz^- 

In the 2D-case, such an operator can be transformed to the the Ginsburg- 
Landau-Zeeman operator (24) by choosing fj, = eB/2hc and multiplying the 
whole operator with —h? /2m. In general dimensions, number k = k/2 means 
the number of particles, and, endomorphisms jz and —jz in the above for- 
mulas are attached to systems electrons resp. positrons. More precisely, by 
the classification of H-type groups, these endomorphisms are acting on the 
irreducible subspaces M"' and the system is interpreted such that there are 
ni/2 particles of the same charge orbiting on complex planes determined by 
the complex structures j^^ resp. — Jz„ in constant magnetic fields whose 
directions are perpendicular to the complex planes where the particles are 
orbiting. The actuality of the complex planes where the orbiting takes place 
can be determined just probabilistically by the probability amplitudes de- 
fined for such systems. The total number of particles is k = {a + b)ni/2. The 
probability amplitudes must refer to k number of particles, that is, they are 
defined on the complex X-space C defined by the complex structure J^^. 
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This theory can be estabhshed just after developing an adequate spectral 
theory. 

Above, adjective "perpendicular" is meant to be just symbolic, for axis Z 
is actually separated from the orbiting. That is, it is not the actual axis of 
orbiting, but it is a vector in the Z-space which is attributed to the orbiting 
by the linear map J : Z — > J^. Thus, it would be more appropriate to 
say that this constant perpendicular magnetic field, B, is just "felt" by the 
particle orbiting on a given complex plane. Anyhow, the B defines a unique 
inertia system on the complex plane in which E = 0, that is, the associated 
electric field vanishes. In all of the other inertia system also a non-zero E 
must be associated with B. The relativistic time T defined on this unique 
inertia system is the inner time defined for the particle orbiting on a complex 
plane. This time can be synchronized, meaning that common time T can 
be introduced which defines time on each complex plane. 

Note that this operator contains also an extra constant term, which is 
explained later as the total energy of neutrinos accompanying the electron- 
positron system. This energy term is neglected in the original Ginsburg- 
Landau-Zeeman Hamiltonian. Thus the higher dimensional mathematical 
model really represents a system of particles and antiparticles which are 
orbiting in constant magnetic fields. Operator, D^», associated with mag- 
netic dipole resp. angular momentum operators, are defined for the lattice 
points separately. Therefore this model can be viewed such that it is asso- 
ciated with magnetic-dipole-moment-crystals, or, angular-moment-crystals. 
In short, they are called Z-crystals. They are particularly interesting on a 
group iJg"'^^ where the Z-space is M.^. On this Euclidean space all possible 
crystals are well known by classifications. It would be interesting to know 
what does this mathematical classification mean from physical point of view 
and if these Z-crystals really exist in nature? 



5.2 Explicit spectra of Ginsburg-Landau-Zeeman operators. 

On Z-crystals, the spectral investigation of the total operator (32) can be 
reduced to the operators induced by A on the invariant subspaces W^. 
On Heisenberg-type groups this operator involves only a single parameter 
^ > 0, where /i^ is the single eigenvalue of — J^. By this reason, it is denoted 
by 

This problem is traced back to an ordinary differential operator acting on 
radial functions, which can be found by seeking the eigenfunctions in the 
form F{X) = /((X,X))h("''^)(X), where / is an even function defined on 
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R and h("''^)(X) is a complex valued homogeneous harmonic polynomial of 
order n, and, simultaneously, it is also an eigenfunction of operator iD^» 
with eigenvalue m/x. Such polynomials can be constructed as follows. 

Consider a complex orthonormal basis, B = . . . , B^}, on the com- 
plex space defined by the complex structure J = (l//x)J^. The corre- 
sponding complex coordinate system is denoted by {zi, . . . , z^}- Functions 

P = zf . . . z^'^zf ...tt satisfying pi H h Pk = P, 91 H V Qk. = n - p 

are re*''-order homogeneous polynomials which are eigenfunctions of \D» 
with eigenvalue m = 2p — n. However, these polynomials are not har- 
monic. In order to get the harmonic eigenfunctions, they must be ex- 
changed for the polynomials 11^^ (P), defined by projections, nl^\ onto the 
space of n*'^-order homogeneous harmonic polynomials of the X-variable. 
By their explicit description (71), these projections are of the form 11^^^ = 

+ sj^'^x^Ax + B^'^^x'^Aj^ + ..., where A^^ = id. By this formula, also 
the harmonic polynomial obtained by this projection is an eigenfunction of 
iD» with the same eigenvalue m/i. 

When operator <l^ is acting on F{X) = f {{X , X))E^^'''^\X) , it defines an 
ordinary differential operator acting on /. Indeed, by Z)^ • / = 0, we have: 



(43) (<^F)(X) = (4(X, X)ri{X, X)) + i2k + 4n)/'((X, X)) 

-(2m/x + 4/.2((l + X))f{{X, X))))H('^'-)(X). 

The eigenvalue problem can, therefore, be reduced to an ordinary differential 
operator. More precisely, we get: 

Theorem 5.1. On a Z- crystal, Br x T\ under a given boundary condition 
Af^B?) + Bf{R^) = defined by constants A,B £ R, the eigenfunctions 
of <iJL can be represented in the form f {{X , X))']\^'^'"^\X) , where the radial 
function f is an eigenfunction of the ordinary differential operator 

(44) (0,,?/)(t) = ^if"{t) + {2k + An) fit) - (2m/i + 4^^2(1 + \i))f{l). 

For fixed degrees n and m, the multiplicity of such an eigenvalue is the 
dimension of space formed by the spherical harmonics H^"'"^) (X) . 

In the non-compact case, when the torus bundle is considered over the 
whole X-space, the L'^ -spectrum can explicitly be computed. Then, the above 
functions are sought in the form f^it) = u{ifi)e~^^^'^ where u{t) is a uniquely 
determined r*"^ -order polynomial computed for fi = 1. In terms of these 
parameters, the elements of the spectrum are i^(f_i^r,n,m) = — ((4r + 4p + A:)/i + 
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This spectrum depends just on p 
element appears for distinct degrees n. 
eigenvalue is infinity. 



= {m + n) /2 and the same spectral 
Therefore, the multiplicity of each 



Proof. Only the last statement is to be established. We proceed, first, with 
the assumption fi = 1. Then, function e~2* is an eigenfunction of this radial 
operator with eigenvalue — (4p + fc + 4). The general eigenfunctions are 
sought in the form 

(45) f{i) = u{i)e-¥. 

Such a function is an eigenfunction of if and only if u[t) is an eigenfunction 
of the differential operator 

(46) {P{^,=i^n,m.)u){i) = Uu"{i) + {2k + 4n - 4t)n'(t) - (4p + A: + 4)n(t). 

Because of differentiability conditions, we impose u'(0) = on the eigen- 
functions. Since in this case, ti(0) ^ hold for any non-zero eigenfunction, 
also the condition u(0) = 1 is imposed. 

In the compact case, corresponding to the ball x torus- type manifolds de- 
fined over a ball -B/j, the spectrum of this Laguerre-type differential operator 
can not be explicitly computed. For a given boundary condition (which can 
be Dirichlet, u{R) = 0, or Neumann, u'{R) = 0) the spectrum consists of 
a real sequence < /xi > > • • • — > — oo. The multiplicity of each of 
these Laguerre-eigenvalues is 1 and the multiplicity corresponding to the 
Ginsburg-Landau-Zeeman operator is the dimension of the space of spheri- 
cal harmonics H^"'*") . The elements of the Laguerre-spectrum are zeros of a 
holomorphic function expressed by an integral formula [CoH]. 

Contrary to the compact case, the spectrum can be explicitly computed 
for the non-compact torus bundle T\H, defined over the whole X-space. An 
elementary argument shows that for any r G N, operator (46) has a uniquely 
determined polynomial eigenfunction 

(47) lt(^=l,r,n,m)(i) = + OiF"^ + 02^"^ H h ttr-lt + 

with coefficients satisfying the recursion formulas 

(48) Oo = 1 , Oi = -ai-i{r - i){r + n + ^k + 1 - i)r~^ . 

Actually, this argument can be avoided and these polynomials can ex- 
plicitly be established by observing that they are nothing but the Laguerre 
polynomials defined as the r'^'-order polynomial eigenfunctions of operator 

(49) Aa{u){i) = tu" + {a + l-t)u', 
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with eigenvalues — r. This statement follows from identity 

(50) Pif.=l,n,m) = 4A( - (4p + /c + 4), 

implying that the eigenfunctions of operators (46) and (49) are the same 
and the eigenvalue corresponding to (47) is 

(51) = -(4r + 4p + /c + 4) , p=^{m + n). 

We also get that, for fixed values of k,n,m (which fix the value also for p), 
functions ii(;i=i,r,n,m) , = 0, 1, . . . oo form a basis in L^([0, oo)). 

In case of a single n, the eigenfunctions are sought in the form 

(52) 7x^™^((X,X))e-|^<^'^>H('^''")(X). 
It turns out that 

(53) «(/i,r,n,m.) (*) = '"(/i=l,r,n,m) (M^) 

and the corresponding eigenvalue is 

(54) i^(/.,r,n,m) = -((4r + 4p + k)fl + 4//2). 

This statement can be explained as follows. For a general fi, the action 
of (55) on a function f(t) = u{^t)e~2^^ can be described in terms of r = /it 
as follows: 

(55) (L(^,„,^)/ = ^{Arfrr + {2k + 4n)/, - {2m + r)/) - 4/^2/, 

from which the statement follows. □ 

This technique extends to general 2-step nilpotent Lie groups, where the 
endomorphisms may have distinct eigenvalues, {/Xi}. In this case the eigen- 
functions are represented as products of functions of the form 

(56) F(,)(X)=/(,)((X,X))H(-"™')(X), 

where the functions in the formula are defined on the maximal eigensubspace 
corresponding to the parameter /ij. When the spectrum is computed on 
the whole X-space, this method works out for the most general Ginsburg- 
Landau-Zeeman operators. In case of a single /i, this method applies also to 
computing the spectra on torus bundles over balls and spheres. In case of 
multiple /i's, it applies to torus bundles over the Cartesian product of balls 
resp. spheres defined on the above Xj-spaces. 
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6 Particles having interior. 

In order to sketch up a clear map for this rather complex section, we start 
with a review of the main mathematical and physical ideas this exposition 
is based on. These ideas are rigorously established in the subsequent sub- 
sections. 

6.1 A preliminary review of the main ideas. 

Systems of particles having insides can be attached to ball x ball- and ballx- 
sphere-type manifolds. Originally they emerged in the second type of spec- 
tral investigations performed in [Sz2]-[Sz4]. These manifolds are defined by 
appropriate smooth fields of Z-balls resp. Z-spheres of radius i?z(x) over 
the points of a fixed X-ball Bx whose radius is denoted by Rx- Note that 
radius Rz{^) depends just on the length, x := |X|, of vector X G Bx over 
which the Z-balls resp. Z-spheres are considered. The centers all of the balls 
resp. spheres which show up in this definition are always at the origin of the 
corresponding spaces. The boundaries of these manifolds are the so called 
spherexball- resp. sphere x sphere- type manifolds, which are trivial Z-ball- 
resp. Z-sphere-bundles defined over fixed X-spheres of radius Rx- In short, 
one considers Z-balls resp. Z-spheres instead of the Z-tori used in the pre- 
vious constructions of Z-crystals. In the isospectrality investigations these 
compact domains corresponding to Rx < oo play the primary interest. In 
physics, however, the non-compact bundles corresponding to Rx = oo (that 
is, which are defined over the whole X-space) become the most important 
cases. In what follows, both the compact and non-compact cases will be 
investigated. 

Contrary to the Z-crystal models, the computations in this case can not 
be reduced to a single endomorphism. Instead, they always have to be estab- 
lished for the complete operator M = ^ daD^* which includes the angular 
momentum endomorphisms Jz with respect to any Z-directions. This oper- 
ator strongly relates both to the 3D angular momentum P = (Pi, P21 -P3) = 
j^Z X p, defined in (21), and the strong interaction term igx'^A'^ta of the 
QCD Lagrangian (25). Actually, it has a rather apparent formal identity 
with Pauli's intrinsic spin Hamiltonian ^ BiPi defined by magnetic fields 
B = {Bi,B2,B3) (cf. [P], Volume 5, pages 152-159). In 5f/(3)-theory the 
term corresponding to the 3D angular momentum is exactly the above men- 
tioned strong interaction term. Due to the new form (26) of the de Broglie 
waves, where the angular-momentum-axes are separated from the planes on 
which the particles are orbiting, also the angular momentum operator has 
to appear in a new form. 
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In order to make the argument about the analogy with Pauh's intrinsic 
spin resp. strong interaction term more clear, note that the Lie algebras 
determined by the 3D angular momenta resp. Gell-Mann's matrices ta are 
su{2) resp. su{3). Thus, the nilpotent-group-models corresponding to these 

classical Yang-Mills models are H^""'^^ resp. the nilpotent group constructed 
by the representations of su{3). When operator M acts on wave function 
(26), it appears behind the integral in the form \Dk*- Consider, first, group 

i/J'^'') and suppose that K = ei, where {61,62,63} is the natural basis on 
M^. Then, on the {62, esj-plane, which is a complex plane regarding the 
complex structure Jg^, operator i-Dei* is nothing but the first component, 
—Pi = i(Z2^^ — Zsg^), of Pauli's angular momentum operator. That 
is, the analogy between M and the classical angular momentum operator 
becomes apparent after letting M act on wave functions (26). This action 
is the very same how P is acting on the original de Broglie waves. Thus this 
new form of action, which can be described by axis-separation and placing 
the orbiting particles onto the complex planes, can really be considered as 
adjustment to the new forms of the wave functions. These arguments work 
out also for groups constructed by su(3)-representations. Thus it is really 
justified to consider M as a spin operator appearing in a new situation. 
The greatest advantage of this new form is that it describes also the strong 
nuclear forces. 

This complication gives rise to a much more complex mathematical and 
physical situation where both the exterior and the interior life of particle 
systems exhibit themself on a full scale. First, let the physical role of the 
Fourier transforms appearing in the formulas be clarified. If term involving 
time is omitted from the formula of wave functions, the rest is called time- 
less probability amplitude. Both these amplitudes and wave functions could 
have been defined also by means of the inverse function e~^^^'^\ Wave func- 
tions (or amplitudes) obtained from the very same function by using e^^^'^^ 
resp. e^'^"^'^^ in their Fourier transforms are said to be wave functions (or 
amplitudes) defined for particle- resp. antiparticle-systems. In other words, 
the definition of particles and antiparticles is possible because of these two 
choices. Calling one object particle and its counter part antiparticle is very 
similar to naming one of the poles of a magnet north-pole and the other 
one south pole. Since the Laplace operators on 2-step nilpotent Lie groups 
are defined by means of constant magnetic fields, this is actually the right 
physical explanation for choosing e^^^'^'^ or e~^^^'^'^ to introduce probabil- 
ity amplitudes. By this reason, M is called unpolarized magnetic dipole 
moment or angular momentum operator. The polarized operators appear 
behind the integral sign of the Fourier integral formula when M is acting 
on the formula. 
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The complexity of this operator is fascinating. For instance, it is the sum 
of extrinsic, L, and intrinsic, S, operators which do not commute with each 
other. Furthermore, operator (E = Ax + (1 + ^x^)A^ + L is a Ginsburg- 
Landau-Zeeman operator which exhibits just orbital spin. The intrinsic life 
of particles is encoded into S. Also the strong nuclear forces, keeping the 
particles having interior together, can be explained by this operator. In 
order to make this complicated situation as clear as possible, we describe, in 
advance, how certain eigenfunctions of A can explicitly be computed. These 
computations provide a great opportunity also for a preliminary review of 
the general eigenfunction computations, which will be connected to these 
particular computations as follows. 

Since these particular functions do not satisfy the boundary conditions, 
they do not provide the final solutions for finding the eigenfunctions yielding 
also given boundary conditions. This conditions can be imposed just after 
certain projections performed on the center (that is, in the insides of the 
particles). But then, these projected functions will not be eigenfunctions 
of the complete A any more. They are eigenfunctions just of the exterior 
operator (E. By this reason, they are called weak force eigenfunctions. The 
strong force eigenfunctions, defined by the eigenfunctions of the complete A 
satisfying also a given boundary condition, can be expressed just by com- 
plicated combinations of the weak force eigenfunctions, meaning that the 
strong forces are piled up by weak forces. Since the weak force eigenfunc- 
tions are Ginsburg-Landau-Zeeman eigenfunctions, by the combinations of 
which the strong force eigenfunctions can be expressed, this theory really 
unifies the electromagnetic, the weak, and the strong nuclear forces. Let 
it be mentioned yet that the only force-category missing from this list is 
the gravitational force. At this early point of the development, we do not 
comment the question if this unification can be extended also to this force. 

Now we turn back to establish an explicit formula describing certain eigen- 
functions of A on general Heisenberg type Lie groups hI^'^'^K More details 
about the general eigenfunction computations will also be provided. Al- 
though this construction can be implemented also on general 2-step nilpo- 
tent Lie groups, this more complicated case is omitted in this paper. First, 
the eigenfunctions of a single angular momentum operator D^*, defined for 
a Z-vector K are described as follows. For a fixed X-vector Q and unit Z- 
vector Ku = ^K^, consider the X-function Qq{X, K^) = {Q + iJ/<„((5),X) 
and its conjugate @q{X, Ku)- For a vector K = k-K^ of length k, these 
functions are eigenfunctions of Dk* with eigenvalue — ki resp. ki. The 
higher order eigenfunctions are of the form QqQq with eigenvalue — p)ki. 
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In order to find the eigenfunctions of the compound operator M^, con- 
sider a Z-sphere bundle ^^^(x) over the X-space whose radius function 
Rz{x) depends just on \X\ = x. For an appropriate function (/)(x, K) (de- 
pending on X and K € Sr^ , furthermore, which makes the following integral 
formula well defined) consider 

(57) TQpgR,mX,Z)= <f e'<^'^>(/.(x,K)(e^e^)(X,i^Jdif„o, 

where dKno is the normalized measure on ^^^(x). By ^ = ^'iDk*, 
this function restricted to the Z-space over an arbitrarily fixed X-vector 
is an eigenfunction of with the real eigenvalue {p — q)Rz{^)- These 
functions are eigenfunctions also of with eigenvalue i?^(x). Also note 
that these eigenvalues do not change by varying Q. 

The function space spanned by functions (57) which are defined by all 
possible (^'s is not invariant with respect to the action of Ax, thus the 
eigenfunctions of the complete operator A do not appear in this form. 
In order to find the common eigenfunctions, the homogeneous but non- 
harmonic polynomials Gq0q of the X- variable should be exchanged for 

the polynomials n^'*(0Q0Q), defined by projections, IIx, onto the space 
oi n = [p + g)-order homogeneous harmonic polynomials of the X-variable. 
Formula IIx = A^-FBi|XpAx + S2x''A^-F. . . , established in (71), implies 
that, over each X-vector, also 

(58) ■HJ'Q„R,mX,Z)= I e'<^'^)0(x,K)nf)(e^e^)(X,if„))dK„, 

Jsr^ 

are eigenfunctions of M and A^ with the same eigenvalues what are defined 
for (57). 

The action of the complete Laplacian is a combination of X-radial differ- 
entiation, 9x, and multiplications with functions depending just on x. Due 
to the normalized measure dKno, these operations can be considered such 
that they directly act inside of the integral sign on function (/>(x, K) in terms 
of the x-variable, only. That is, the action is completely reduced to X-radial 
functions and the eigenfunctions of A can be found in the form 

(59) /(x') / e'<^'^>FQp,(X,K„))dK„o, where 

JSr^ 

(60) Fqp,{X,K^)) = ^{K)U^;^\%{X,Ku)Ql{X,Ku)). 

The same computations developed for the Z-crystals yield that this reduced 
operator appears in the following form: 

(61) (0^(i),t/)(t) = 4t7"(t) + (2A: + 4n)/'(t)-(2m/z(t) + V(t)(i + it))/(t), 
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where t = x^, and = Rz{Vt) = Rzi'x.). Note that this is exactly the 
same operator what was obtained for the radial Ginsburg-Landau-Zeeman 
operator on Z-crystals. Since function fj,{i) may depend also on i = x^, 
it actually appears in a more general form here. However, for constant 
radius functions Rz, it becomes the very same operator, indeed. That is, 
also this eigenfunction-problem is reduced to finding the eigenfunctions of 
the Ginsburg-Landau-Zeeman operator reduced to X-radial functions. This 
reduced operator remains the same by varying Q, thus also the spectrum on 
the invariant spaces considered for fixed Q's is not changing regarding these 
variations. This phenomena reveals the later discussed spectral isotropy 
yielded on these models. 

Note that this construction is carried out by a fixed X- vector Q, but it 
extends to general polynomials as follows. Consider an orthonormal system 
B = {Bi, . . . , -Br} of vectors on the X-space. They form a complex, but 
generically non-orthonormal basis regarding the complex structures Jku> 
where the unit vectors yielding this property form an everywhere dense 
open set on the unit Z-sphere. This set is the complement of a set of 
measure. The corresponding complex coordinate systems on the X-space 
are denoted by {zk^i = ®Bk^^,---,zk^^, = Qbk^J- For given values 
pi, gi, . . . consider the polynomial HiLi ^^Kui^KuV Then functions 

(62) I e'(^'^^f{^'MK)uP n ^'K^^^A.idK^o = 

-'^Rz 1=1 

= /(x2) / e'<^'^)FBp,g.(X, Ku)dKno = n^Bp^g^BzimX, Z) 

are eigenfunctions of A if and only if function /(x^) = f[i) is an eigenfunc- 
tion of the radial operator (61), where p = "^PijQ = "^qijn = p + q. 

Consider a Z-ball bundle with radius function Rz (x) defining a compact or 
non-compact ball x ball- type domain. Then the eigenfunctions satisfying the 
Dirichlet or Z-Neumann conditions on this domain can not be sought among 
the above eigenfunctions because functions FQpq{X, K^) resp. F^piqi {X, Ky) 
are not spherical harmonics regarding the i^T^-variable but they are rather 
combinations of several spherical harmonics belonging to different eigenval- 
ues of the Z-spherical Laplacian. It turns out, however, that one can con- 
struct the complete function space satisfying a given boundary condition 
by the above formulas if functions F,, are substituted by their projections 
n^^(F..) into the space of s^'^-order spherical harmonics regarding variable 

Ku. Actually, this projection appears in a more subtle form, H^"*^ = H^^, 
which projects i^,,,, first, into the space of (i' + a)*''-order homogeneous poly- 
nomials (where v resp. a refer to the degrees of functions to which '^{K) 
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resp. the {p + q)*''-order polynomials are projected). The projection to 
the s*''-order function space applies, then, to these homogeneous functions. 
The most important mathematical tool applied in this investigations is the 
Hankel transform developed later. 

Although these new functions yield the boundary conditions, they do 
not remain eigenfunctions of the complete Laplacian anymore. However, 
they are still eigenfunctions of the partial operator CE. They are called 
weak-force-eigenfunctions, which can be considered as electromagnetic-force- 
eigenfunctions because both ffi and the Ginsburg-Landau-Zeeman operators 
can be reduced to the same radial operator. The old ones from which the new 
functions are derived are called linkage-eigenfunctions of A, which bridge the 
electromagnetic interactions with the weak interactions. 

Since the extrinsic, L, and intrinsic, S, operators do not commute, the 
weak-force-eigenfunctions can not be the eigenfunctions of the complete 
operator A = (E + S. In other words, the weak-force-eigenfunctions can 
not be equal to the strong-force-eigenfunctions defined by the eigenfunc- 
tions of A satisfying a given boundary condition. To construct these func- 
tions, the Fourier integrals must be considered on the whole Z-space M', 
by seeking them in the form f^i e'^^'^^(^Q(x, z)n^^(F,,(X, i^))(iEr. The ac- 
tion of A on these functions can be described in the form J^^i e'^^'^'* Qo 

((/>!, . . . , (^rf)n^^ K))dK, where operator Oa(0ii • • • > 4'd), correspond- 

ing d-tuples of (X, Z)-radial functions to each other, is defined in terms of 
Hankel transforms combined with radial derivatives of the functions appear- 
ing in the arguments. Finding the eigenfunctions of A means finding the 
eigen-d-tuples, ((/>!, . . . , cpd), of the radial, so called roulette operator Qa- As 
it will be pointed out, these eigenfunctions exhibit properties characteristic 
to strong force eigenfunctions. 

These arguments really unify the 3 fundamental forces of particle theory. 
The real union is exhibited by the common unpolarized operator A. After 
polarization, they are separated into three categories. These cases corre- 
spond to the function spaces on which the polarized operators are acting. 
The details are as follows. 

6.2 Twisted Z-Fourier transforms. 

This is the main mathematical tool which incorporates de Broglie's wave 
theory into the new models in a novel, more general form. The name in- 
dicates that the Fourier transform is performed, over each X-vector, only 
on the Z-space in the same manner as if one would like to consider the de 
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Broglie waves only in the center of the Lie group. However, an important 
new feature is that this transform apphes to product of functions, where 
one of them purely depends just on the center variable, K, while the other 
is a complex polynomial of the X-variable defined in terms of the complex 
structures Jru^ where = A'/k and k = \K\. This Z-Fourier transform 
is said to be twisted by the latter polynomials. Thus this transform has 
impact also on the X-variable. This simple idea establishes the necessary 
connection between the abstract axes, Ku, and the particles placed onto the 
complex planes of the complex structures Jk^ ■ 

This transform is defined in several alternative forms corresponding to 
those introduced in the previous review. The difference is that, over each 
X-vector, the following functions and integrals are defined on the whole Z- 
space. This is contrary to the previous section where the integral is defined 
just on Z-spheres. Since the eigenfunctions constructed in the review do not 
satisfy any of the boundary conditions, this reformulation of the Z-Fourier 
transform is really necessary for the complete solutions of the considered 
problems. 

In the first case, consider a fixed X-vector Q, and define the same functions 

(63) eQ{x,K^) = {Q + ijKAQ),x),eQ{x,K^), 

as above. For fixed integers p,q >0 and L^-function (j){x,K), consider the 
Z-Fourier transform 

(64) :FQpqmX, Z)= [ e'<^'^>(/.(x, K)e^(X, Ku)el{X, Ku)dK, 

which is said to be twisted by the polynomial Qq{X, Ku)QQi-^j Ku)- Func- 
tion (j) is considered also in the form (/)(x, k)(^(A„), where (/9 is a homogeneous 
polynomial of the K-variable. The L^-space spanned by these functions is 
denoted by ^g^g) where n = p+q indicates that these functions are ?i*''-order 
polynomials regarding the X-variable. The space spanned by the twisted 
functions (pix, K)ePQ{X, Ku)Qq{X, Ku) is denoted by P^Qpg. This is the 
pre-space to which the Fourier transform is applied. 

Instead of a single vector Q, the second alternative form is defined regard- 
ing K = k/2 independent vectors, B = {Ei, . . . , of the X-space. Such a 
system forms a complex basis for almost all complex structure Jr-^, where 
Ku = K/]i. Now the twisting functions are polynomials of the complex 
coordinate functions 

(65) {zkAX) = eQ,iX,Ku), . . .,ZK^^{X) = eQ^ix,Ku)}, 

where these formulas indicate that how the coordinate functions can be 
expressed in terms of the above ©g-functions. For appropriate functions 
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(/)(x, K) and polynomial exponents {pi,qi) (where i = 1, . . . ,k/2) transform 
^Bipm^'P^i^^^) is defined by: 



where ip{Ku) is the restriction of an m*^-order homogeneous polynomial, 
ip{K), onto the unit sphere of the K-space. If X^j(pi + Qi) = n, then the 
twisting functions are n*''-order complex valued polynomials regarding the 
X- variable and for any fixed X, the whole function is of class regard- 
ing the K-variable. These properties are inherited also for the transformed 
functions. 

When cl){'x.,K) runs through all functions which are of class L^, for any 

(n) 

fixed \X\, the transformed functions span the function space ^p>„ „ . All 
these function spaces, defined for index sets satisfying n = '^{pi + qi), span 
the function space denoted by = X]{(p^g.)} *bp,i2j- '^^^ corresponding 
pre-spaces are denoted by P^gp g, and P$g^ respectively. 

Twisted Z-Fourier transforms (57)-(62) defined in the previous section 
by considering Dirac-type functions concentrated on spheres Sr of radius 
-R(x) can be generated by the familiar L°° approximation of these Dirac 
type functions by functions. That is, this radius depends just on x. 
Function spaces ^t^I^w ^bw ^^'^ their pre-spaces P*Bp,gi_R ^^sp. P$g]j 
are defined similarly as for functions. They are defined also for one-pole 
functions. These cases are denoted such that B is replaced by Q. 

In many cases the theorems hold true for each version of these func- 

(n) 

tion spaces. By this reason we introduce the unified denotation $ ^ and 

(n) 

P$ jj, where the dots represent the symbols introduced above on the indi- 
cated places. A unified denotation for the total spaces are ^^'^ and F^^^ , 
indicating that symbols piQi do not show up in these formulas. If letter R 
is omitted, the formulas concern the previous cases when the functions are 
of class regarding the if-variable. 

The elements of the latter total function spaces are complex valued func- 
tions defined on the (X, Z)-space such that, upto multiplicative X-radial 
functions, they are n*''-order polynomials regarding the X- variable and for 
any fixed X they are functions regarding the Z-variable. It is very im- 
portant to clarify the relations between the above twisted spaces and the 
latter complex valued function spaces. By considering an arbitrary real 
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basis Q = {Qi, ■ ■ ■ , Qk} on the X-space, the complex valued functions ap- 
pear in the form E{ai,...,afc} '^ai,...,afc (|^|, 2') nf=i(Qi> ^)''S where the sum 
is considered for all sets {ai, . . . ,ak} of non-negative integers satisfying 
J^Oj = n and functions (p are L^-functions for any fixed \X\. The func- 
tion space spanned by these functions is called straight L\ space of complex 
valued (X , Z) -functions. The Z-Fourier transform performed on such func- 
tions is called straight Z-Fourier transform. The explanations below show 
that twisted functions (pi\X\,K) Y\i=i Zx^ii^P^^ii^) 

can be converted to 

straightly represented functions in the terms of which the twisted Z-Fourier 
transform becomes a straight Z-Fourier transform. This conversion works 
out also in the opposite (from the straight to the twisted) direction. The 
precise details below show that the straightly defined function spaces are 
complete regarding the norm in which both $g and the pre-space 
are everywhere dense subspaces. By this reason, the corresponding straight 
spaces are denoted by resp. P^b- However, these spaces are equal. Ul- 
timately, both versions of the Z-Fourier transforms define authomorphisms 
(one to one and onto maps) of this very same ambient function space. 

This situation can be illuminated by the real k x k matrix field Aij{Ku) 
defined on the unit Z-vectors which transforms the real basis Q to the vector 
system Br = {Si, . . . , 5^/2, ^(fc/2)+i = Jku{Bi), . . . , Bk = J/c^ (-8^/2)}- 
That is, this field is uniquely determined by the formula Bi = Yl^=i ^ijQji 
where i = 1, . . . , fc. The entries are polynomials of K^. By plugging these 
formulas into the twisted Z-Fourier transform formula, one gets the straight 
representations both of the twisted functions and their Z-Fourier transforms. 

Conversion from the straight to the twisted functions is more complicated. 
In this case vectors Qi should be exchanged for vectors Bj according to the 

formula Qi = '}2^j=i ^ij^^j- Then vectors Bj resp. 5(fc/2)+j; where j < k/2, 
should be expressed in the form 

(67) Bj = ^{B, + JKAB,)) + ^iB,-JKABj)) resp. 

(68) B(fc/2)+,- = JkABj) = \{B, + JkSBj)) - \{B, - JkABj)), 

which, after performing powering and appropriate rearranging, provide the 
desired twisted formulas. Due to the degeneracy of the matrix field Aij on 5b, 
some entries of the inverse matrix field A~/ have limits -|-oo or —00 of order 
at most k/2 on this singularity set. Therefore, the Z-Fourier transform of 
a term of the twisted function involving such functions may be not defined, 
despite the fact that the Z-Fourier transform of the whole function exist 
which is equal to the transform of the straightly represented function. In 
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other words, the infinities appearing in the separate terms cancel each other 
out in the complete function. 

This contradictory situation can be resolved as follows. For a given e > 0, 
let Sse be the e-neighborhood of the singularity set on the unit Z-sphere 
and let MSbc be the conic set covered by the rays emanating from the origin 
which are spanned by unit Z- vectors pointing to the points of S-Be- For 
an L|-function K) discussed above let (^^(x, be the function which 
is the same as (j) on the outside of MS'se and it is equal to zero in the 
inside of this set. Then, regarding the L^-norm, lim^^QCp^ = cp holds. For 
a function F{X,K), expressed straightly, define F^{X,K) by substituting 
each function (p by (p^. If function tp^{K) is defined by 1 outside of M^Be 
and by inside, then F,(X, K) = tp^{K)F{X, K) holds. Convert into the 
twisted form. Then the twisted Z-Fourier transform is well defined for each 
term of the twisted expression, providing the same transformed function as 
what is defined by the straight Z-Fourier transform. Thus each straightly 
represented function is an L^-limit of functions which can be converted 
to twisted functions in which each twisted term is an function having 
well defined twisted Z-Fourier transform. In this process, function ip^(K), 
which is constant in radial directions, can be chosen such that it is of class 
satisfying < tpeiK) < 1, furthermore, it vanishes on M5'b^/2 ^^'^ is 
equal to 1 outside of MSbe- Thus we have: 



Theorem 6.1. The twisted functions, defined in each of their terms by 
functions, form an everywhere dense subspace P$g in the complete space 
P$B straightly defined functions. Although the first space depends on 
B, the second one is a uniquely determined function space which does not 
depend neither on B nor on Q. 

A general function from the ambient space becomes an appropriate twisted 
function belonging to the dense subspace after multiplying it with a function 
ipe{K) which is zero on the above described set MSbe CLnd equal to 1 on the 
complement of this set. This function, ipe{K), can be chosen to be of class 
C°° such that it is constant in radial directions satisfying < ipe{K) < 1, 
furthermore, it vanishes on ^i^d is equal to 1 outside of KS'bc ■ Then 

the L^- approximation is defined by the limiting e ^ 0. 

The twisted Z-Fourier transform continuously extends to the straight Z- 
Fourier transform defined on the ambient space. On the ambient space this 
Fourier transform is an authomorphism. The Z-Fourier transform is a bi- 
section between the everywhere dense subspaces and Particularly 
the relation P$b = holds. 
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The same statements hold true also for Q-pole function spaces. Total 
spaces resp. P$g^ are spanned by the corresponding Q-pole functions 
whose poles, Q, are in the real span of vectors belonging to B. 

In order to show the one to one property, suppose that the Z-Fourier 
transform of a twisted function vanishes. That is: 

k /2 

(69) / e^<^'^^> Yl '^fe,..)(l^l'^)ri-l.(^)^l.(^)rf^ = 

{{P,,q.)} 




Then, for any fixed X, the Z-Fourier transform of g^)} ^{pt,qi){^j 

vanishes. Therefore, this function must be zero for all X and for almost 
all K. For vectors Ku not lying in the singularity set Sb, the B is a com- 
plex basis and complex polynomials Yl^^i z^^-{X)'Zj^^-{X) of the X- variable 
are linearly independent for distinct set {{pi,qi)} of exponents. Thus all 
component functions <j)(^p.^q.^{\X\, K) of the K variable must vanish almost 
everywhere. This proves that the twisted Z-Fourier transform is a one to 
one map of the everywhere dense twisted subspace onto the everywhere 
dense twisted range space. The same proof, applied to straight functions, 
yields the statement on the straight space where the Z-Fourier transform is 
obviously an onto map by the well known theorem of Fourier transforms. 

This section is concluded by introducing several invariant functions by 
which physical objects such as charge and volume will be precisely defined. 
Suppose that system B consists of orthonormal vectors and let B"*- be the 
orthogonal complement of the real span, B = Span^(B), of vectors belong- 
ing to B. Basis B defines an orientation on B. Let B^ resp. B:;*: be the 
two possible orientations which can be chosen on B"*" such that, together 
with the orientation of B, they define positive resp. negative orientation 
on the X-space, M^'. Charging a particle system represented by B means 
choosing one of these two orientations. Once the system is charged, choose 
an orthonormal basis also in B"*-, complying with the chosen orientation. 
Together with B, this system defines an orthonormal basis, Q, on the X- 
space by which the matrix field Aij{Zu) defined above can be introduced. 
The invariants of these matrix field are independent from the above ba- 
sis chosen on B"*-. Invariant functions ch{Zu) = Tr{Aij{Zu)) — {k/2) and 
v{Zu) = det{Aij{Zu)) are called charger and volumer respectively. Their in- 
tegral regarding a proper probability density defines the charge resp. mass 
of the particle being on a given proper state (these concepts are precisely 
described later in this paper). 
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6.3 Hankel transform. 

This twisted Fourier transform is investigated by means of the Hankel trans- 
form. The statement regarding this transform asserts: 

Theorem 6.2. The Fourier transform considered on M' transforms a prod- 
uct, f{r)F^'^\9), of radial functions and spherical harmonics to the product, 
Hv\f){r)F'^'^\0) , of the same form, i. e., for any fixed degree v of the 

spherical harmonics, it induces maps, Hu\f){r), on the radial functions, 
which,so called Hankel transform, is uniquely determined for any fixed in- 
dices I and V. 

This is actually a weak form of the original Hankel theorem which can 
be directly settled by the following mean value theorem of the spherical 
harmonics, F^'^\9), defined on the unit sphere, S*, about the origin of M'. 
These functions are eigenfunctions of the Laplacian with eigenvalue X^, 
moreover, there exists a uniquely determined radial eigenfunction (fx^{p), 
where < p < vr and fx^{0) = 1, on S which has the same eigenvalue 
Xi, such that, on a hypersphere crp{6) C S of radius p and center 9 on the 
ambient sphere S, the identity j^F^^^duno = F^''\e)ipx,{p) holds, where 
dcFno is the normalized measure measuring o" by 1. 

This mean value theorem can be used for computing the Fourier transform 
j^i e'^^^^^f{\K\)F(''\eK)dK at a point Z = {\Z\,ez). This integral is com- 
puted by Fubini's theorem such that one considers the line Iz^t), spanned 
by 9z', parameterized with arc-length t; and satisfying lz{Z) > 0, and one 
computes the integrals first in hyperplanes intersecting Iz at t perpendic- 
ularly and then on Iz by dt. On the hyperplanes, write up the integral 
in polar coordinates defined around the intersection point with Iz, where 
the radial Euclidean distance from this origin is denoted by r. Consider 
polar coordinates also on S around 6z, where the radial spherical distance 
from this origin is denoted by p. Then r = |t||tan/9| holds, furthermore, a 
straightforward computation yields: 

(70) [ e'^^^''^ f{\K\)F^^\eK)dK = Hl^\f){r)F^^\ez), where 

e"-'\tt' / f(\ttanp\)^xApf^^^dpdt, 

-oo Jo I COS PI 

where il/_2 denotes the volume of an (I — 2)-dimensional Euclidean unit 
sphere. These formulas prove the above statement completely. 
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6.4 Projecting to spherical harmonics. 

Among the other mathematical tools by which twisted Z-Fourier transforms 
are investigated are the projections Il^^'^\ip{Ku) Yli=i co'^" 
responding s*''-order polynomials to r*^-order polynomials of the ETu-variable. 
Although the functions they are applied to may depend also on the X- 
variable, they refer strictly to the -fCu- variable, meaning, that they are per- 
formed, over each X-vector, in the Z-space. These characteristics are exhib- 
ited also by the fact that these projections appear as certain polynomials of 
the Laplacian Ak^ defined on the unit K-sphere. It is much more convenient 
to describe them in terms of homogeneous functions, which are projected 
by them to harmonic homogeneous polynomials. In this version these poly- 
nomials depend on K and the projections can be described in terms of the 
Laplacian Ak defined on the ambient space. By restrictions onto the unit 
K-spheres, one can easily find then the desired formulas in terms of n^'*\ 

For an n*'^ order homogeneous polynomial, Pn{K), projection 11^^ : = 

nj^'"^ onto the space of n*^ order harmonic polynomials can be computed 
by the formula 

(71) uPiP^iK)) = Y,Ci^HK,KrA^j,{P^{K)), 

s 

(n) 

where Cq = 1 and the other coefficients can be determined by the recur- 
sive formula 2s(2(s + n) — l)ci"'* -|- Cg"\ = 0. In fact, exactly for these 
coefficients is the function defined by an arbitrary Pn{K) on the right side 
a homogeneous harmonic polynomial. These formulas can be easily estab- 
lished for polynomials Pn{K) = (W,i^)", defined by a fixed Z-vector W. 
Since they span the space of n'^^ order homogeneous polynomials, the state- 
ment follows also for general complex valued polynomials. This projection 
is a surjective map of the ?i*^-order homogeneous polynomial space 
onto the space, 7i^'^\K), of the n*''-order homogeneous harmonic polyno- 
mials whose kernel is formed by polynomials of the form {K , K) Pn-2{K) ■ 
Subspace H^") C is a complement to this kernel. 

The complete decomposition of Pn{K) appears in the form Pn{K) = 
Y^^{K, KYHPn-2i{K), where i starts with and running through the inte- 
ger part, [?T./2], of n/2, furthermore, HPn-2i{K) are harmonic polynomials 
of order {n — 2i). It can be established by successive application of the 
above computations. In the second step in this process, one considers the 
functions Pn{K) — n^^(P„(K)), which appear in the form {K, K)Pn-2{K), 
and obtains HPn~2{K) by the projection 11^" '^\Pn^2{K)). Then, also 
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this second term is removed from Pn{K), in order to get ready for the 
third step, where the same computations are repeated. This process can 
be completed in at most [n/2] steps. It is clear that projections H^'" 
resulting functions HPn-2i{K) from Pn{K) are of the form 11^'"' = 

D{n,n-2i)^^K ^^'"^ ^'^^X' where indices (n, n — 2i) indicates that the projec- 
tion maps n'^^-order polynomials to (n — 2i)*^-order harmonic polynomials 
(in this respect, projection 11^^ is the same as n^'""^). The technical cal- 
culation of constants -D(n,n-2i) is omitted. The corresponding projections 

n^''*^ defined on the unit spheres can immediately be established by these 
projections defined for homogeneous functions. They appear as polynomials 
of the Laplacian IS.Ku defined on the unit K^-sphere. 

Since these projections depend just on the degrees r and s, they apply 
also to twisted functions which depend on the X-variable as well. In order 
to make the Hankel transform applicable, they are used for decomposing 
functions in the form 

(72) mi'J{^{K^) n ^Kj^r^tj^)) 

(r;s) «=1 

= /„(x,k)n^jF(P"'?')(x,i^„)), 

where the right side is just a short way to describe the sum appearing on 
the left side in terms of compound indices a = (r, s) and functions F^^''''-\ 

To be more precise, functions (/'n(x, K)P^'^\X, K^) appearing in the pre- 
spaces must be brought to appropriate forms before these projections can 
directly be applied to them. First of all, function </)„(x, K) should be con- 
sidered in the form (j)n{'x.,K) = (/)„_t,(x, k)^?^'') (J^), where Lp^^\K) is 
an ?;*'*-order homogeneous harmonic polynomial. Then, after implement- 
ing all term- by-term multiplications in the products of = {Bi,X) + 
{iJKu{Bi),X) and = {Bi,X) - {{Jx^iBi), X) , the above polynomials 
have to be taken to the form p(")(X,K„) = ELo ^^"'"^(^' -^^«)' where 
polynomial involves exactly a number of linear polynomials of the 

form (Jx^iQi), X) . The above projections defined in terms of r directly act 
on functions Lp^'"\K)P^"''°'\X , K) satisfying r = v + a, which are obtained 
by term by term multiplications of the sums given above for (j)n{x,K) and 
p(^\X,Ku). 

This complicated process can be considerably simplified by considering 
only one-pole functions defined for single Q's which are in the real span of 
the vector system B. In fact, all the 1-pole total spaces = Yip q^Qpg^ 
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where n = p + q, span also the total space g.) *b^' ^^^^ there is really 
enough to establish the theorems just for these simpler one-pole functions. In 
this case, function P^"''"'\X, Ku) is nothing but a constant-times of function 

(73) i?j''^)(x,ir„) = {Q,xr-'^{jKAQ),xr = {Q,xr-'^{[Q,xiKur. 

Note that, depending on p and g, the component of a particular p(")(X, Ku) 
corresponding to a given < a < n may vanish. However, there exist such 
pairs {p,q) for which this a-component is non-zero. Space P("'°)(X, IT) 
is defined by the span of functions Rq''°'\x,K). Note that has been 
changed to K in the above formula. These functions are n^^- resp. a*^-order 
homogeneous polynomials regarding the X- resp. K- variables. Keep in mind 
that these functions can be derived from functions QqQq, where p + q = n, 
by linear combinations, therefore, they belong to the above twisted function 
spaces. More precisely, for a given n, there exist an invertible matrix M^'°^ 
such that r'^q'''\x, K^) = Y.p,q ^S'^^^q^q hold, where p + q = n. 

Projections = d"\^, ^)'A5f acting on functions P{X,Ku) = 
%{X,Ku)Wq{X,Ku) resp. P{X,Ku) = Hti where 
q = n resp. ^iPi + <?« = n, regarding the X- variable are also involved 
to these investigations. Since they are n^'^-order homogeneous functions 

(n) 

regarding X, projection 11^ applies to them immediately. Then, for any 

fixed Ku, function n5^^(P(X, iiT^)) is an 71 -order homogeneous harmonic 
polynomial regarding the X-variable. The twisted Z-Fourier transforms, 
TCJ^s(^p.q.^{(p){X, Z), involving these projections are defined by 



(74) / e'^^^^U{\X\,K)uP([[z^^jX)z%jX))dK. 

The corresponding function spaces spanned by the transformed func- 
tions and the pre-space are denoted by 3^^^^, and PSBpig, respectively. 
These function spaces are well defined also for one-pole functions and also 
for the third type of Z-Fourier transforms defined for Dirac type generalized 
functions concentrated on Z-sphere bundles. 

Although the following results are not used in the rest part of this paper, 
because of their importance, we describe some mathematical processes by 
which these projections can explicitly be computed. Further on, the for- 
mulas concern a fixed n even if it is not indicated there. Decomposition 
into K-harmonic polynomials will be implemented for the K-homogeneous 
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functions (p^'"\K)Rq {X , K) of order r = v + a. According to two different 
representations of the first function, these projections will be described in 
two different ways. The first description is more or less technical, yet, very 
useful in proving the independence theorems stated below. In the second 
description, the projected functions are directly constructed. In both cases 
we consider Q-pole functions defined by a unit vector Q G Spank's = B. 
However, the multipole cases referring to vector systems B are also discussed 
in the theorem established below. 

According to the formula (74), a pure harmonic one-pole function with 
pole ( in the K-space is of the form if^^\K) = D2sk^^(C) -f^)''"^'') where 
-^0 = iCr s-rid the other coefficients depending on ICT"^* and constants Cs 
are uniquely determined by the harmonicity assumption. It is well known 
that these functions span the space of homogeneous harmonic polynomials. 
Functions d'^^'^^ obtained by directional derivatives regarding a fixed vector 
K are also homogeneous harmonic C-pole functions of order (v — c). 

The action of operator A^- on ip^^'\K)RQ\x, K) results function: 



c=0 

The terms of this sum are obtained such that acts on Rq''"'\x, K), 
resulting the very last term, while the others are due to the action of A''"'^ 
on the product according to the formula d^'^^j^^^ d^~'^RQ\x, K). Note that, 

because of the harmonicity, the action of A^ on (J"^^ is trivial. 

is) 

When the complete projection 11^'^ is computed, then the 6's involved to 
the formula are denoted by hj. For given c, factor out from the 

corresponding terms. Thus the final projection formula appears in the form 
^jK|^=Pi'^(X,i^)|[Q,X]|2c, where term Pi'^ is equal to 

(75) ^Z)e,(Q,X)""''(a^'^-V['^))(ir)(Jc(Q),X)^^-^(J^(Q),X)'^-^^-^ 

3 

This is a rather formal description of the projected functions. A more 
concrete construction is as follows. The linear map X ^ Z defined by 
X — > is surjective whose kernel is a (fc — /)-dimensional subspace. 
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Next projections Hk will be investigated in the Z-space over such an X- 
vector, X, which is not in this kernel and the unit vector (^q(X) defined 
by = |[(5,X]|Cq(X) is not in the singularity set 5b- If Xq denotes 

the orthogonal projection of X onto the Z-dimensional subspace spanned by 
Q and vectors Jx^iQ) considered for all unit vectors Ku, then = 

The direct representation of (^(^') is the product of two harmonic one-pole 
functions having perpendicular poles. One of the poles is C = Cq{^) while 
the other is an arbitrary perpendicular unit vector (^-^ 5b- Then, for all 
< c < V, consider the product ^p^^^ = ip^^^ip^y of harmonic homogeneous 

one-pole functions- Because of the perpendicular poles, these products are 
also harmonic functions of the K- variable, furthermore, considering them for 
all < c < v and C"*", they span the whole space of v^'^-order homogeneous 
harmonic polynomials of the K-variable, for each point X. However, is not 
pointing into the direction of [Q, X] in general. When such a pure K-function 

is multiplied with Rq'"'\x, K), over X, for the required projections, only the 

decomposition of ip^~' "^^C, A')" should be determined. A simple calculation 

shows: 

(76) ^'^y-''\c,Kr = j2D,{K,KY^^y-''^''-''\ 

i=0 

This function multiplied with {Q, X)"'~"'\[Q, X]\°'ip^y^ provides the desired 
decomposition and projections. When this projected function is considered 
over an arbitrary X, neither (p^^l nor ip^^ appear as pure K-functions. 

One can state only that, for all i, this product is an a^'^-order homogeneous 
polynomial, also regarding the X-variable. This complication is due to that 
Cq{X) and Cq{X) are not parallel, thus the projection operator involves 
both functions to the computations. According to these arguments we have: 

Theorem 6.3. (A) For given non-negative integers r, a, and {v — a) < 

s < {v -\- a); where s has the same parity as {v + a) or {v — a); func- 
tions <p^^\K)R^q''^\x,K) project to {y -\- a)*^ -order K-homogeneous func- 
tions which, over X , appear in the form 



(77) 

where < c < f , which, by omitting term {K, K)^^^"""^^^"^ , are s^^-order 
homogeneous K-harmonic polynomials. For given v and a, this projection is 
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trivial (that is, it maps to zero) for all those values s which do not satisfy 
the above conditions. Above, exactly the non-trivial terms are determined. 

(B) Functions {K, K)^'"~^"'~'^^^'^ and of degrees (v + a — s) resp. 
(n — a) are independent from c, furthermore, function ip^^^ip^- is an a^^- 
order homogeneous polynomial of the X-variable. It follows that, for a given 
s, the subspaces, PHoq''^\ spanned by projected functions considered for dis- 
tinct pairs {v, a) are independent, more precisely, ^^.^ PHog'"'* = PHog^ 
is a finite direct sum decomposition of the corresponding space of complex 
valued homogeneous functions which are n^^-order regarding the X-variable 
and s^^-order harmonic functions regarding the Z-variable. This decompo- 
sition is further graded by the subspaces, PHog''^''^'' defined for distinct c 's 
introduced above. 

( C) These statements remain true for the total spaces PHo^'"^ and PHo^ 
obtained by summing up all the corresponding previous spaces defined for Q 's 
which are in the real span of independent vector-system B. Since the projec- 
tions Hx and Hk commute, they remain true also for spaces PXog'"^ and 

PXOg^ (as well as for versions defined for fixed Q's) obtained by applying 
Hx to the corresponding spaces PHo. 



Proof. The independence-statement in (B) follows also from (75), because 
functions [[QjXjp'^ regarding distinct powers 2c are independent and, for a 
fixed Ku, also functions Pc'^\x,K) of the X-variable expressed by means 
of complex structure as a complex valued function has distinct real 
and imaginary degrees with respect to distinct a's. That is, these functions 
defined for a fixed a can not be a linear combination of the others defined 
for other a's. 

Statement (C) can be established by an appropriate generalization of (75). 
Instead of Q-pole functions, now functions belonging to PHog'"'^ should be 
projected, where a = ^ aj and degree regards Bi. In this situation, one 
gets functions {[Bi,X],[Bj,X])'"'^ multiplied with the corresponding func- 
tions Pc^j {X,V), which, for a fixed and system Cij of exponents have 
distinct real and imaginary degrees regarding with respect to distinct 

Ei^i = a 7^ a' = X^iO^- □ 
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6.5 Twisted Hankel decomposition. 

In order to construct the complete pre-spaces P^q\ P$g\ P^q\ and 

P$b\ by means of functions in PHog'*''"^ resp. PHog'*''"\ they must 
be multiphed with functions of tlie form </>(|X|, \K\) which multiplied with 
|^|i)+a incorporated to the PHo-functions provide K-radial L|^-functions 
for any fixed \X\ = x. These functions summed up regarding (f,a) span 
the corresponding n^'^-order twisted spaces whose L^-closures provide the 
complete space which can also be introduced by straightly defined functions. 
Since they are everywhere dense in the straightly defined function spaces, 
the twisted spaces must be complete regarding the PHo-spaces. Thus the 
completion of twisted spaces is ultimately implemented on the space of K- 
radial functions. 

Projections 11^^ should be defined for functions belonging to p^^''^'") 

resp. p^^''''")^ first. Note that this operation has no effect on radial func- 
tions. It's action is restricted to the spherical harmonics defined by restrict- 
ing the above homogeneous K-harmonic polynomials to the unit K-sphere 
over each point X. By this interpretation, these projection operators can be 
expressed as polynomials of the Laplacian A/^^ defined on this unit sphere. 
The function spaces obtained by projecting the whole corresponding ambi- 
ent spaces are denoted by P^q'^^ resp. P$g'*\ There is described in the 
previous theorem that which functions labeled by {v, a) provide non-trivial 
s-components in these operations. They provide direct sum decompositions 
of the ambient spaces which is called also pre-Hankel decomposition. When 
also (which operator commutes with for all s) is acting, the ob- 
tained pre-Hankel spaces are denoted by PSq'^^ resp. PS^^'^\ where n still 
indicates the degree of the involved homogeneous polynomials regarding the 
X-variable. 

The Z-Fourier transforms / e'<^'^>P*J^''^di^ resp. / e'^^'^^P^^'^^dK 

define the s^'^ -order twisted Hankel spaces H$q''*^ resp. H^g''*^ which are 

projected to HHg'*^ resp. HHg'**^ by The direct sums of these non- 

complete subspaces define the corresponding total twisted Hankel spaces. 
They are different from the pre-spaces, but also everywhere dense subspaces 
of the corresponding straightly defined function spaces. Thus the closure of 
these twisted spaces provides again the whole straightly defined spaces. 
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6.6 Twisted Dirichlet and Z-Neumann functions. 

All above constructions are implemented by using the whole center M'. In 
this section twisted functions satisfying the Dirichlet or Z-Neumann condi- 
tion on the boundary, dM, of ball x ball- type domains, M, are explicitly 
constructed by the method described in the review of this section. That is, 
they are represented by twisted Z-Fourier transforms of Dirac type general- 
ized functions concentrated on dM. Due to this representation, the eigen- 
functions of the exterior operator (E satisfying given boundary conditions 
can explicitly be computed. 

In the first step of this process consider a sphere, Sr, of radius R around 
the origin of the Euclidean space R'. As it is well known, the Dirichlet 
or Neumann eigenfunctions of the Euclidean Laplacian — Ajgz on the ball 
Bji bounded by Sr appear as products of s*^-order spherical harmonics 
if^^^Ku) with radial functions y^'^\z). For s = 0, these eigenfunctions are 
radial taking 1 at the origin and having multiplicity 1. For s > 0, the radial 
functions take at the origin and the multiplicity, for fixed s and i, is equal 
to the dimension of the space of s^'^-order spherical harmonics (p^^\Ku)- 
Corresponding to the Dirichlet or Neumann conditions, these eigenvalues 

(s) (s) 

are denoted by A]-,- and X)^- respectively. For any fixed s and condition D 
or A^, these infinite sequences satisfy A-*^ t oo and, except for = X^^\, also 

(s) 

the relation < A- holds. 

(s) 

Eigenfunctions corresponding to Dirichlet or Neumann eigenvalues A,- = 
A can be represented by the integral formula 



(78) yf^(z)(^(^)(Zj= / e'<^'^)(^(^)(K/k)(iK, 



where dKno is the normalized integral density on the sphere of radius \/A- 
Apply — A]g! on the right side to see that this function is an eigenfunction 

(s) 

of this operator with eigenvalue A- := A, which, because of the Hankel 
transform, must appear as the function being on the left side. This formula 
strongly relates to the third version of the twisted Z-Fourier transforms 
which is introduced also in the review of this section. However, it can 
be directly applied there only after decomposing the functions behind the 
integral sign by the projections II^'''^ 
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Before this application, functions y^'^\z) := y{t), belonging to an eigen- 
value A, are more explicitly determined as follows. By formulas 

(79) Az = dtdt + ^-^dt + ^, As^^'^ = -s{s + l-2)ip(^\ 

it satisfies the differential equation 

// ^ - 1 / s(s + l- 2)^ 

(80) y" + -^y' + {X- -^—^^ '-}y = 0, 

which, after the substitutions r = \/At and y{t) = z{t), becomes 

„ I — 1 , r s(s + I — 2) ^ 
T 

That is, function J(t) = T'/^~^z(r) satisfies the ordinary differential equa- 
tion 

(82) J" + h' + {1 - l^i±l^} J = 0, 

therefore, it is a bounded Bessel function of order [s + 1/2 — 1). Thus, except 
for complex multiplicative constant, equation J = Js+«/2-i niust hold. 

In order to find the functions satisfying the boundary conditions, consider 
spheres 5';^(x2) of radius A(x^) around the origin of the Z-space. For appro- 
priate functions </>(x^) and ip{Ku), the twisted Z-Fourier transform on the 
sphere bundle Sx is defined by: 

(83) TQp,x{<t)^){X,Z)= i S^^'^^(t>{i^)^{Ku){%%){X,Ku)dKno. 

where dKno = dK/Vol[ax) is the normalized measure on the sphere and 
function ip{K /\^)QPqWq{X , K /\^) is defined on Sx- Since no Hankel pro- 
jections are involved, these functions do not satisfy the required boundary 
conditions, yet. However, by the above arguments we have: 

Theorem 6.4. Consider a ballxball-type domain defined by the Z-balls 
-S_R(x2) cif^d let \^^\^^) := A(x^) be a smooth function defined by the i^^ 
eigenvalue in the s^^ -order Dirichlet or Neumann spectra of the Euclidean 
balls Bji(^2y Then function (83) defined for A = A^*''(x^), or 

(84) / e'<^'^),/)(|X|2)vrW(/3Mn5^)(G^e^Q))(X,KjdK„, 

satisfy the Dirichlet resp. Z-Neumann condition on the domain M . The 
restrictions of these functions onto M span $q'*^(M) resp. 's!'q^\m), for 
any fixed boundary condition. If functions ©qQq «re exchanged for the 
polynomials HiLi ii-^)'^K A^), then the above construction provides the 
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Dirichlet-, Z-Neumann-, resp. mixed- condition- functions spanning 
resp. 3^"'^)(M). 

6.7 Constructing the orbital and inner force operators. 

The complicated action of the compound angular momentum operator 
on twisted Hankel functions is due to fact that its Fourier transform, iDx*, 
is non-commuting with the Hankel projections, that is, the commutator on 
the right side of Dx • = Jlf^^Dx • +[Dk • n^^] is non- vanishing in 
general. Thus there are two non-trivial terms on the right side of equation 

(85) Mz I e'<^'^>/a(x,k)n^jF(P-'^-)(X,i^„))dif = 
= i J e'^^'''^>U^,k)DK •Ilf,^{F'j'-'>^\X,Ku))dK = 

= -^j e^<^'^>/«(x,k)(n^^Dx.+[D^.,n^j)(F(^'-''>)(x,i^„))di^, 

by which the orbital: 

(86) Lz j e'<^'^>/,(x,k)n?^jF(f-''»)(X,K„))dK = 

= i j e'^^'""^ U^,mK^DK • {F^^-'^^Hx,K^))dK 
and the intrinsic spin operators: 

(87) Sz j e'<^'^>/,(x,k)n^jF(f-'?')(X,if„))dif = 

= \j e^<^'^^>/,(x,k)[Z)A-,n^J(F(?'-«»)(X,K„))dJf 
are defined, respectively. 

The commutator appears in the following explicit form 

(88) [DK^.IiU = S'^K^Mj,^ = M;,x5^n^^, 

where Mj^± = Mk—Q^Dk^* is the Z-spherical angular momentum operator 
defined on Z-spheres. (For a fixed unit vector Ku and orthonormal basis 
{ea,Ku} (where a = l,2,...,/ — l)of K-vectors, this operator is of the form 
M^x = "^^daDa*. ) Formula (88) follows by applying relations Dk • 
= Ak^Dk • — 2M^x and the commutativity of Ax„ with operators 
da,dK^, Da», and Dk^* in order to evaluate 

n . Tt('''''~2*) — n - n •n{r-2i,r-2i) 
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where, according to the denotations introduced above, r = v + a holds. This 
computation resuhs the equation [Dk»,II^^'^ ^*^] = Pj+i(Ax„)Mj^x, where 
the lowest exponent of Aj<-^ in the polynomial Pj_|_i(Ai^^) is i + 1. This term 
defines a uniquely determined constant times of projection 2(i+i)) g^^j^ 

that P,+i{AkJ = ^i+in5^'J"^^*+^^^ + Pi+2i^Ku) holds, where the lowest 
exponent of A.Ku in the polynomial Pi^2{^Ku) is i + 2. In the next step, the 
above arguments are repeated regarding to obtain Pj+3. This process 
can be concluded in finite many steps which establish the above formula 
completely. 

Formula (88) allows to define an inner algorithm where these computa- 
tions are iterated infinitely many times. In the second step it is repeated 
for fj^^ := S'^fa as follows. First note that operator /^S'^II^^M^x = 
M^±n^^5'^/a can be decomposed in the following form: 

(89) n^^M^x = -dKMf)K.PK. • +M^n^jf . 

The action of the first operator on functions 

(90) F^P'''\x,Ku) = ip^[\K^)@l{X,K^)el{X,K^) 

k/2 

(91) resp. F(f'«)(X,i^„)=v.(^-)(A'.)n4„.^i„„ 

i=l 

where ip^^^K^^) denotes an r*^-order homogeneous harmonic polynomial in- 
troduced at explicit description of Hankel projections, results 

This term is already in finalized form which does not alter during further 
computations. Together with the orbiting operator, they define, in terms of 
the radial operator 

(92) 0^aHfp^,■ • • , /&) = -{P - q)ma + ^i^JZ/^^f )), 

the one-turn operator 

(93) M« I e^<^'^)/,(x,k)n^jF(^'-5')(X,K„))dir = 

= / e'<^'^^ Oi'^ {U,,...,U,WK^{F^P-''^\x,K^))dK, 

where the name indicates that it is expressed in terms of the first power of 
the inner spin operator Sa permutating the Hankel radial functions. 

Be aware of the novelty of this spin-concept emerging in these formulas! 
It defines just permutation of the radial Hankel functions to which no actual 
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spinning of the particles can be corresponded. This concept certainly does 
not lead to a dead-end-theory like those pursued in classical quantum theory 
where one tried to explain the inner spin of electron by actual spinning. 

This abstract merry-go-round does not stop after making one turn. It 
is actually the second operator, Mi^-II^ fj^\ (^9) which generates the 
indicated process where the above arguments are repeated for functions 

(2) 

:= Spfa- Index 2 indicates that these functions are obtained from the 
starting functions fa^ '■= fa by the next step. The details are as follows. 

In these computations the second operator is derived by the Hankel trans- 
form turning functions defined on the -fsT-space to functions defined on the 
K-space. Operator denotes the inverse of the Hankel transform Hg^ , 

where denotes the third index in /? = (v, a, s). Then we have: 



where fg = '^'^HIb )■ this step commutator [Dj>»,Il^ ] can be 

(3) ~(2) 

calculated by (88), resulting functions := S^fa which are subjected to 
the operations performed in the following step 3. These steps must infinitely 
many times be iterated. 

~(2) 

Regarding radial functions fa , which are defined on the K-space, the 
orbiting spin and the one-turn operator can be defined in the same way 
as they are defined on the K-space. After performing the Z-Fourier and 
the associated Hankel transforms on these finalized functions they become 
functions defined on the iT-space. By adding these terms to those obtained 

(2) 

in the first step, one defines the two-turn operator associated with 

(2) 

Qa {fpii ■ ■ ■ 1 ff}^), called two-turn merry-go-round and roulette operators 
respectively. The sum of the orbiting spin operators defines L^^^, which is 
called one-turn orbiting spin operator involving just S'^ into its definition. 
One should keep in mind that these operators involve also the operators 
defined in the previous step. These computations work out for an arbitrary 

it'^^-step defining operators Q)^a\fpi,---,fi3ci)^ ^'^'^ L*^") which are 

called U-turn merry-go-round-, roulette-, and {u — l)-turn orbiting spin op- 
erator respectively. In the end, the action of = can be described 
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in the form: 

(95) Mz I e'<^'^>/,(x,k)n^jF(^'"^')(X,i^,))di^ = 

where operator Oa(//3i> • • • > //3 J = Oa°°^(//3i, • ■■Jf3a), called high roulette 
operator, is defined by the infinite series lim^^oo OcT^ • In this sense, the 
complete angular momentum operator can be called high merry-go- 
round operator. The point in this formula is that this action can be described 
in terms of d-tuples, (//j^ , . . . , /^^), of radial functions which can not be 
reduced to a single function defined by a fixed a. 

This statement holds also for the total Laplacian A, thus the eigenfunc- 
tions satisfying the Dirichlet or Z-Neumann conditions on the considered 
manifolds can also be completely described in terms of radial functions. 
But the corresponding equations involve all radial functions (//3i , ■ ■ ■ , 
defined for all indices By this reason, operator Qo pl^ys the role of a 
confider, giving rise to the effect that the eigenfunctions satisfying a given 
boundary condition are expressed in terms of radial functions which do not 
satisfy these conditions individually but exhibit them together, confined in 
a complicated combination. Neither can these eigenfunctions be observed as 
single ffi-force potential functions. On the other hand, for a given bound- 
ary condition, the (E-eigenfunctions span the whole L^-space, therefore, A- 
eigenfunctions satisfying the same boundary condition can be expressed as 
infinite linear combinations of CE-eigenfunctions. This observation further 
supports the idea that the rather strong A-forces are built up by the much 
weaker ffi-forces in a way how Hawking describes the action of the Weinberg- 
Salam weak force in [H], pages 79: 

"The Weinberg-Salam theory known as spontaneous symmetry breaking. 
This means that what appear to be a number of completely different particles 
at low energies are in fact found to be all the same type of particle, only in 
different states. At high energies all these particles behave similarly. The 
effect is rather like the behavior of a roulette ball on a roulette wheel. At 
high energies (when the wheel is spun quickly) the ball behaves in essentially 
only one way - it rolls round and round. But as the wheel slows, the energy 
of the ball decreases, and eventually the ball drops into one of the thirty- 
seven slots in the wheel. In other words, at low energies there are thirty 
seven different states in which the ball can exist. If, for some reason, we 
could only observe the ball at low energies, we would then think that there 
were thirty-seven different types of ball!" 
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This quotation explains the origin of the name given to the roulette opera- 
tors. This polarized operators are derived from the non-polarized merry-go- 
round operators, which name was suggested to me by Weinberg's book [Wl] 
where the name "merry-go-round" appears in a different situation not dis- 
cussed here. All these arguments clearly suggest that the A-eigenfunctions 
must correspond to the strong forces keeping the quarks together. A for- 
mula expressing these eigenfunctions as linear combinations of weak force 
potential functions would shed light to the magnitude of the strong force 
piled by the roulette operator during building up these eigenfunctions. 

The explicit eigenfunction computations give rise to difficult mathemati- 
cal problems which are completely open at this point of the developments. 
In this paper we explicitly describe only the eigenfunctions of the exterior 
orbiting operator ffi defined by omitting the interior spin operator from the 
A. This operator is a scalar operator whose action can be reduced to a sin- 
gle radial function /. By the arguments developed in the introductions and 
also at several points in the body of this paper, this operator is associated 
with the weak force interaction. These computations are established in the 
next section. This section is concluded by saying more about the recovery 
of several concepts of the standard model within this new theory. 

Suppose that function 



is an eigenfunction of the complete Laplacian A. It is called also probability 
amplitude of the particle system. Also in this formula the Einstein conven- 
tion indicates summation. A fixed index a = {v, a, s) is called slot-index 
and function 



is the so called high slot probability amplitude. These exact mathemat- 
ical objects correspond to quarks whose fiavor is associated with index v 
and its color is associated with index a. Index s is called azimuthal index. 
Note that also these indices have mathematical meanings, namely, they re- 
fer to the degrees of the corresponding polynomials by which the formula 
of a slot-amplitude is built up. According to these definitions, a particle- 
system-amplitude is the sum of the high slot probability amplitudes which 
are considered to be the mathematical manifestations of quarks. Slot am- 
plitudes defined in the same way by the CE-eigenfunctions are called retired 
or laid-off slot probability amplitudes. 

It is very interesting to see how these constituents of a high energy par- 
ticle are held together by the strong force interaction. When Dk* is acting 
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on functions {Q, X) resp. {Jk{Q)i ^) of a quark, then the first one becomes 
of the second type and the second one becomes of the first type. In either 
cases the color index, a, changes by 1 and a slot-ampUtude of odd color 
index becomes an amplitude of even color index. This process can be inter- 
preted as gluon exchange in the following way. Action of Dk* on {Q,X) 
resp. {Jk{Q), X) are interpreted as gluon absorption resp. emission. More 
precisely, some of the slot-particles (quarks) of odd color index emit gluon 
which is absorbed by some of the slot-particles (quarks) which also have 
odd color index. As a result they become quarks of even color index. The 
same process yields also for quarks having even color index, which, after 
gluon exchange, become quarks of odd color index. This is a rather clear 
explanation for the flavor-blindness and color sensitiveness of gluons. 

It also explains why a high slot-particle (quark) can never retire to be- 
come a laid-off (E-particle. Indeed, in a high slot-state defined for fixed slot 
index a = {v, a, s) the corresponding Hankel function does not satisfy the 
boundary conditions, however, it can be expanded by the (E-eigenfunctions. 
At this point nothing is known about the number of ffi-eigenfunctions by 
which a high slot-amplitude can be expressed. This number can very well 
be equal to the infinity, but it is always greater then one. Then, instead of 
consisting of a single term, the high slot probability density is a multi-term 
sum of weak densities determined by these laid off probability amplitudes. 

Real positive function ipip whose integral on the whole ball x ball- type do- 
main is 1 is called probability density. Protons are defined by functions 
ch{Zu){ipc'4'c)i^^ where ■0c is a constant times of ip, whose integral on 
the ball X ball- type domain is 1. If this integral is or —1, it is called neu- 
tron or antiproton respectively. (The mass can be defined by means of 
\det{Aij{Zu))\, but we do not go into these details here.) This argument 
shows that in a complete eigensubspace decomposition of the function 
space of a two-step nilpotent Lie group representing a particle system the 
eigenfunctions, actually, represent all kind of particles and not just particu- 
lar ones. This phenomenon can be considered as a clear manifestation of the 
bootstrap principle, from which super string theory grew out, also in the new 
theory. In super string theory the notion was (cf. [G], pages 128) that "a set 
of elementary particles could be treated as if composed in a self-consistent 
manner of combinations of those same particles. All the particles would 
serve as constituents, all the particles (even the fermions in a certain sense) 
would serve as quanta for force fields binding the constituents together, and 
all the particles would appear as bound states of the constituents" . 
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6.8 ffi-forces in the union of 3 fundamental forces. 

^(n) 

The eigenfunction computations of CE on the twisted function space 3 
satisfying a given boundary condition can be reduced to the same radial 
differential operator what was obtained for the standard Ginsburg-Landau- 
Zeeman operator on Z-crystals. To see this statement, let the complete 
Laplacian (32) act on (84). By the commutativity relation <f = jiDx* 
and IS.Z § = — ^ k^, this action is a combination of X-radial differentiation, 
9x, and multiplications with functions depending just on x, that is, it is 
completely reduced to X-radial functions. More precisely, 

(98) /^{'HrQp.Rm) = <f e^<^'^>0ij,x2(<^)/3nJ)(e^e^)(iK„o 

J Sr 

holds, where 

(99) (0ii,x2</')(x2) = 4x2</,"(x2) + {2k + 4(p + g))</.'(x2) + 

-(i22(l + ix2) + (p-g)i?)0(x2), 

and (j)' , (f)" mean the corresponding derivatives of (t){t) with respect to the t 
variable. The eigenfunctions of A can be found by seeking the eigenfunctions 
of the reduced operator 0r,x2 among the X-radial functions. 

(s) 

Note that no projections 11]^ are applied in the above integral formula, 
thus these eigenfunctions do not satisfy the boundary conditions in general. 
These projections, however, do not commute with Dk*, and the eigenfunc- 
tions of the complete operator A can not be expressed in terms of a single 
function </>(x^). This simple reduction applies just to the exterior opera- 
tor, CE, defined by neglecting the anomalous intrinsic momentum operator 
from A and keeping just the orbital (alias, exterior) spin operator L^. 
Regarding this operator we have: 

Theorem 6.5. The exterior operator ffi, on constant radius Z-ball bundles 
reduces to a radial operator appearing in terms of the Dirichlet-, Neumann-, 
resp. mixed-condition-eigenvalues \\^^ of the Z-ball Bz{R) in the form 

(100) (0,,-/)(t) = 4t7"(t) + (2fe+4n)/'(t)-(2my^+4^(l + -t))/(t). 

By the substitution X = \J A^*^/4, this is exactly the radial Ginsburg-Landau- 
Zeeman operator (55) obtained on Z-crystal models. 

Despite this formal identity with electromagnetic forces, the nuclear forces 
represented by ffi manifest themself quite differently. Like for the electro- 
magnetic forces, one can introduce both charged and neutral particles also 
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regarding ffi. A major difference between the two particle- systems is that the 
particles represented by (E are extended ones. This property can he seen, for 
instance, by the eigenf unctions of A which involve also Z-spherical harmon- 
ics. Due to these harmonics, the multiplicities of the eigenvalues regarding 
ffi are higher than what is corresponded to the same eigenvalue regarding 
the Ginsburg-Landau-Zeeman operator on Z-crystals. An other consequence 
of the extension is that the (E-neutrinos always have positive mass, which is 
zero for neutrinos associated with electromagnetic forces. 

Actually, the nuclear forces represented by ffi are weaker than the elec- 
tromagnetic forces. This phenomenon can be explained, by the extension, by 
the very same argument of classical electrodynamics asserting that the elec- 
tromagnetic self-mass for a surface distribution of charge with radius a is 
e^/Q'Kac^, which, therefore, blows up for a — > 0. Note that, in the history of 
quantum theory, this argument provided the first warning that a point elec- 
tron will have infinite electromagnetic self-mass. It is well known that this 
problem appeared with even much greater severity in the problem of infinities 
invading quantum field theory. The tool by which these severe problems had 
been handled is renormalization, which turned QED into a renormalizable 
theory. The above argument clearly suggests that also the (E-theory must be 
renormalizable. 

The major difference between the ffi- and the complete /^.-theory is that 
the first one is a scalar theory, which can be reduced to a radial operator 
acting on a single radial function, while the reduced operator obtained in 
the A-theory acts on d-tuples of radial functions. As it is described above, 
this action defines also a new type of nuclear inner spin of the extended 
particles to which new particles such as quarks and gluons can be associated 
and by which strong nuclear forces keeping the parts of the nucleus together 
can he introduced. Partial operator Qi is the maximal scalar operator in A. 
Except for ffi = (E^^^ , partial operators ffi^") ( defined by replacing by 

L^^J and A^") (defined by replacing by M^^j can be reduced just to 
operators which irreducibly act on d-tuples of radial functions. 

The main unifying principle among the three fundamental forces is that 
they are derived from the very same Hamilton operator, A. More precisely, 
the Hamilton operators of the individual elementary particles emerge on cor- 
responding invariant subspaces of A such that it is restricted to the sub- 
space corresponded to the given elementary particles. The corresponding 
forces are defined by those acting among these particles. The electromag- 
netic forces manifest themself on function spaces consisting functions which 
are periodic regarding a Z-lattice Tz. The systems of particles defined on 
these function spaces are without interior. They consist of particles such as 
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electrons, positrons, and electron-positron-neutrinos. The various nuclear 
forces appear on function spaces defined on Z-ball and Z-sphere bundles by 
fixed boundary conditions. The attached particles, which do have interior, 
and the forces acting among them are discussed above. 

The function spaces corresponded to particular particles are constructed 
with the corresponding twisted Z-Fourier transforms. Since the eigenfunc- 
tions of the Hamilton operators can be sought in this form, this transform 
seems to be the only natural tool for assigning the invariant function spaces 
corresponding to elementary particles. It is also remarkable that the twisted 
Z-Fourier transforms emerge also in the natural generalization of the de 
Broglie waves fitting the new theory. 

A particle-system defined by a fixed function space consists of all kind 
of particles which can be defined by the characteristic property of the given 
function space. For instance, all particles without interior appear on a Tz- 
periodic function spaces. In order to avoid the annihilation of particles by 
antiparticles, the particles belonging to the same system are considered to be 
not interacting with each other. According to this argument, the particles in 
a system defined by an invariant subspace are gregarious which are always 
accompanied with other particles. For instance, an electron is always par- 
tying with an electron-neutrino. This complexity of the particle- systems is 
reflected by the Laplacian (Hamiltonian) which appears as the sum of Hamil- 
tonians of particles partying in a system. This phenomenon is a relative of 
those described by the bootstrap principle of super string theory. 

Let it be mentioned yet that the distinct function spaces belonging to dis- 
tinct type of forces are not independent, thus there is the possibility to work 
out an interaction theory between particles of distinct types. The existence 
of such a viable theory is a completely open question in this field. 

The CE-forces are very similar to the weak nuclear forces described in the 
Weinberg-Salam theory. Weinberg introduces these forces on pages 116-120 
of his popular book [Wl] as follows: 

" The weak nuclear force first turned up in the discovery of radioactivity 
by Henri Becquerel in 1896. In the 1930's it become understood that in the 
particular kind of radioactivity that was discovered by Becquerel , known 
as beta decay, the weak nuclear force causes a neutron inside the nucleus 
to turn into a proton, at the same time creating an electron and another 
particle known today as antineutrino, and spitting them out of the nucleus. 
This is something that is not allowed to happen through any other kind 
of force. The strong nuclear force that holds the protons and neutrons 
together in the nucleus and the electromagnetic force that tries to push the 
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protons in the nucleus apart cannot change the identities of these particles, 
and the gravitational force certainly does not do anything of the sort, so 
the observation of neutrons changing into protons or protons into neutrons 
provided evidence of a new kind of force in the nature. As its name implies, 
the weak nuclear force is weaker than the electromagnetic or the strong 
nuclear forces. This is shown for instance by the fact that nuclear beta 
decay is so slow; the fastest nuclear beta decays take on the average about 
a hundredth of a second; languorously slow compared with the typical time 
scale of processes caused by the strong nuclear force, which is roughly a 
millionth millionth millionth millionth of a second. 

In 1933 Enrico Fermi took the first significant step toward a theory of this 
new force. ... There followed a quarter century of experimental afford aimed 
at tying up the loose ends of the Fermi theory. ... In 1957 this [problem] 
was settled and the Fermi theory of the weak nuclear force was put into its 
final form. ... Nevertheless, even though we had a theory that was capable 
of accounting for everything that was known experimentally about the weak 
force, physicists in general found the theory highly unsatisfactory.... The 
things that were wrong with the Fermi theory were not experimental but 
theoretical. ... when the theory was applied to more exotic processes it 
gave nonsensical results.. ..when they did the calculations the answer would 
turn out to be infinite.... Infinities like these had been encountered in the 
theory of electromagnetic forces by Oppenheimer and others in the early 
1930's, but in the late 1940's theorists had found that all these infinities in 
quantum electrodynamics would cancel when the mass and electric charge 
of the electron are properly defined, or "renormalized" . As more and more 
became known about the weak forces it became increasingly clear that the 
infinities in Fermi's theory of the weak forces would not cancel in this way; 
the theory was not renormalizable. The other thing that was wrong with the 
theory of weak forces was that it has a large number of arbitrary elements.... 

I had worked on the theory of weak forces off and on since graduate school, 
but in 1967 I was working instead on the strong nuclear forces, the forces 
that hold neutrons and protons together inside atomic nuclei. I was trying 
to develop a theory of the strong nuclear forces based on an analogy with 
quantum electrodynamics. I thought that the difference between the strong 
nuclear forces and electromagnetism might be explained by a phenomenon 
known as broken symmetry, which I explain later. It did not work. I found 
myself developing a theory that did not look like at all the strong forces 
as they were known to us experimentally. Then it suddenly occurred to 
me that these ideas, although they had turned out to be completely useless 
as far as the strong forces were concerned, provided a mathematical basis 
for a theory of weak nuclear forces that would do anything that one might 
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want. I could see the possibility of a theory of the weak force analogous 
to quantum electrodynamics. Just as the electromagnetic force between 
distant charged particles is caused by the exchange of photons, a weak force 
would not act all at once at a single point in space (as in the Fermi theory) 
but it would be caused by the exchange of photonlike particles between 
particles at different positions. These new photonlike particles could not 
be massless like the photon (for one thing, if massless they would have 
been discovered long before), but they were introduced into the theory in 
a way that was so similar to the way that the photon appears in quantum 
electrodynamics that I thought that the theory might be renormalizable in 
the same sense as quantum electrodynamics-that is, that the infinities in the 
theory could be canceled by a redefinition of the masses and other quantities 
in the theory. Furthermore, the theory would be highly constrained by its 
underlying principles and would thus avoid a large part of arbitrariness of 
previous theories. 

I worked out a particular concrete realization of this theory, that is, a 
particular set of equations that govern the way the particles interacted and 
that would have the Fermi theory as a low energy approximation. I found 
in doing this, although it had not been my idea at all to start with, that 
it turned out to be a theory not only of the weak forces, based on an anal- 
ogy with electromagnetism; it turned out to be a unified theory of the weak 
and electromagnetic forces that showed that they were both just different as- 
pects of what subsequently became called an electroweak force. The photon, 
the fundamental particle whose emission and absorption causes electromag- 
netic forces, was joint in a tight-knit family group with the other photonlike 
particles predicted by the theory: electrically charged W particles whose ex- 
change produces the weak force of beta radioactivity, and a neutral particle 
I called the "Z", about which more later. (W particles were an old story in 
speculations about the weak forces; the W stands for "weak". I picked the 
letter Z for their new sibling because the particle has zero electric charge 
and also because Z is the last letter of the alphabet, and I hoped that this 
would be the last member of the family.) Essentially the same theory was 
worked out independently in 1968 by the Pakistani physicist Abdus Salam, 
working in Trieste.... Both Salam and I had stated our opinion that this 
theory would eliminate the problem of infinities in the weak forces. But we 
were not clever enough to prove this. In 1971 I received a preprint from a 
young graduate student at the University of Utrecht named Gerard 't Hooft, 
in which he claimed to show that this theory actually had solved the prob- 
lem of the infinities: the infinities in calculations of observable quantities 
would in fact all cancel just just as in quantum electrodynamics...." 
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7 Unified wave mechanics. 

There are two natural ways to furnish the time on nilpotent groups. One of 
them is defined by solvable extensions while in the other case the time-axis is 
introduced by Cartesian product with the real line M. These two extensions 
are called expanding- and static-models respectively. 

The metric on the nilpotent group is positive definite where the Lapla- 
cian turns out to be natural physical Hamilton operator corresponding to 
elementary particle systems. Concrete systems are represented by the corre- 
sponding invariant subspaces of the Laplacian and one obtains the Hamilton 
operator of a given system by restricting the Laplacian onto these subspaces. 
In order to establish the wave equations regarding these Hamiltonians, on 
both extensions indefinite metrics must be introduced. That is, adequate 
wave mechanics associated with these Hamiltonians can just relativistically 
be introduced such that one assumes appropriate Lorenz-indefinite metrics 
on both extensions. As it turns out, these metrics really provide the familiar 
wave operators of wave mechanics. 

7.1 Solvable extensions. 

Any 2-step nilpotent Lie group, N, extends to a solvable group, SN, which 
is defined on the half-space J\f x ]R-|_ by the group multiplication 

(101) {X, Z, t){X',Z', t') = [X + t^X', Z + tZ' + ^t^ [X, X'],tt'). 

This formula provides the multiplication also on the nilpotent group N, 
which appears as a subgroup on the level set t = 1. 

The Lie algebra of this solvable group is S = J\f (BT. The Lie bracket is 
completely determined by the formulas 

(102) [dt,X] = ^X ; [dt,Z]=Z ; [J^,J^]/sN = W,^f]/N, 
where X £ X and Z £ Z. 

The indefinite metric tensor is defined by the left-invariant extension of 
the indefinite inner product, ( , ), defined on the solvable Lie algebra S by 
{dt,dt) = -1, {dt,M) = 0, and = {M,M)m. The last formula 

indicates that the original innerproduct is kept on the subalgebra M. Lie 
algebra S is considered as the tangent space at the origin (0, 0, 1) G SN 
of the solvable group and ( , ) is extended to a left-invariant metric g onto 
SN by the group multiplication described above. Scaled inner product ( , )q 



Z. I. SZABO 



71 



with scahng factor g > is also defined by {dt,dt)q = — but, keeping 
both the inner product on M and the relation dt -L M on M intact. That 
is, the scaling regards just the direction regarding dt- The left invariant 
extension of these inner products are denoted by g^. 

For precise explanations we need some explicit formulas on these groups 
as well. The left-invariant extensions, Yj, Vq,T, of the unit vectors 

(103) Ei = di , ea = da , e = qdt 
picked up at the origin (0, 0, 1) are the vector fields: 

(104) Yi=t5X, ; Va = tZ^ ; T = qtdt, 

where Xj and Zq, are the invariant vector fields introduced on N previously. 

One can establish these formulas by the following standard computations. 
Consider the vectors di,da and dt at the origin (0,0, 1) such that they are 
the tangent vectors of the curves Cyi(s) = (0, 0, 1) + sOa, where A = i,a, t. 
Then transform these curves to an arbitrary point by the left multiplications. 
Then the tangent of the transformed curve gives the desired left invariant 
vector at an arbitrary point. 

According to these formulas, not t but T defined by 5^ = T is the correct 
physical time parameterization on the t- parameter line, which, by the below 
arguments, are geodesies on SN . The transformation law dx = {dt/dT)dt 
yields the relations (dt/dT) = qt, Int = qT, and t = e'^'^ . Thus a t-level set 
is the same as the T = (In t)/g- level set and subgroup N corresponds both 
io t = 1 and T = 0. The reversed time —T is denoted by r. 

Let Cx{s) and Cz{s) be integral curves of finite length ||c^|| resp. \\cz\ \ of the 
invariant vector fields X and Z on the subgroup N. Then the flow generated 
by dr moves these curves to c!^(s) resp. cj(s) of length ||c!JJ|| = 1 1 1 1 e'^'^/^ 
resp. Ilc^ll = llc^lle'^'^. That is, by considering them as functions of the time- 
variable r, the length is increasing such that the rate of change (derivative 
with respect to r) is proportional to the length of the curves. In other words, 
this mathematical space-time model represents an expanding micro universe 
where the distance between particles is growing exactly in the same way how 
this growing distance was measured by Hubble, in 1929, between galaxies 
[H]. 

Edwin Hubble came to this conclusion after experimenting red-shift in 
the spectra of galaxies he observed for cataloguing their distances from the 
earth. It was quite a surprise that, contrary to the expectation, the red- 
and blue-shifted galaxies occurred with no equal likelihood but most of them 
appeared red-shifted. That is, most of them are moving away from us. Even 



72 



EXACT MATH-MODELS OF MICRO UNIVERSES 



more surprising was to find that even the size of a galaxy's red shift is not 
random, but it is directly proportional to the galaxy's distance from us. This 
means, that the farther the galaxy is, the faster it is moving away. This is 
the familiar Hubble's law which was actually predicted by the Friedmann 
cosmological model, in 1922. 

The tendency to expand must be rooted from the very same tendency 
inbuilt into the microscopic universe. This argument, however, contradicts 
the experimental fact (explained in the introduction) according to which 
the particles are not expanding. This paradox can be resolved, however, by 
recognizing that, due to the expansion, the change of the constant magnetic 
fields, which are present in the particles also by the new model, induces 
electromagnetic waves which are immediately radiated out into the space. 
This explanation clarifies not only this paradox but it casts also a new light 
to the presence of the constant radiation experimentally known in the space. 
These arguments greatly enhances the importance of the exact mathematical 
models introduced in this paper. 

The covariant derivative can be computed by the well known formula 

(105) {VpQ,R) = ^{{P,[R,Q]) + {Q,[R,P]) + {[P,Q],R)}, 
where P, Q, R are invariant vector fields. Then we get 

Vx+z{X* + Z*) = V^+ziX* + Z*) - q{^{X, X*) + {Z, Z*))T; 
VxT = |x ; VzT = qZ ; VtX = VtZ = VtT = 0, 

where denotes covariant derivative and X, X* ^ X; Z, Z* ^ Z;T ^ T. 

The Laplacian on these solvable groups can be established by the same 
computation performed on A^. Then we get 

(106) A = t^Az - q^t^dl + tdt) + 

1 ' ' k 

t{Ax + (MX), Jp{X))dlp + Y.d^D^.)+ q\- + l)tdt = 

q:;/3=1 «=1 

= ^'(^Az-d\^^ 
1 ' ' k 

a;/3=l a=l 

This is the Laplacian on the solvable extension of a general 2-step nilpotent 
Lie group. On the extension of a H-type group it appears in the following 
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simpler form: 

(107) A = {e2«^Az-92^} + {e'?^(Ax + ^x2Az+j;a«I)«.)+g(|+/)aT} 

This operator is expressed regarding the cohapsing time direction. Substitu- 
tion T = —T provides the operator in terms of the expanding time direction. 
Note that the first operator, e~^^'^A^ — d^^, looks like "expanding meson 
operator", while the second one is similar to the Schrdinger operator of 
charged particles. This question is further investigated in the next section. 

In order to understand the deeper connections to general relativity, also 
the Riemannian curvature should explicitly be computed. This calculation 
can straightforwardly be implemented by substituting formulas (106) into 
the standard formula of the Riemannian curvature. Then we get 

(108) Rg{X* AX) = R{X* A X) + ^[X* , X] AT + ^X* A X; 

2 

(109) Rq{X AZ) = R{X AZ) + | Jz(X) A T + A Z; 

(110) Rg{Z* AZ) = R{Z* AZ) + q^Z* AZ; 

(111) Rg{{X + Z),T){.)=qV^^^z{.); RgHX + Z) AT) = 

(112) = ig(^ UX) A - Jl) - q^i^X + Z) AT, 

a 

where is the 2-vector dual to the 2-form {Jz {Xi),X2) and R is the 
Riemann curvature on N. The vectors in these formulas are elements of the 
Lie algebra. 

By introducing H{X,X*,Z,Z*) := {Jz{X), Jz*{X*)) , for the Ricci cur- 
vature we have 

(113) Ri^^x) = Ri{X) - q\^ + l^)X- 

(114) Rig{Z) = Ri{Z) - q\^ + l)Z ; Rig{T) = +q\-^ + l)T, 
where the Ricci tensor Ri on is described by formulas 



Ri{X,X*) = -i^ii'(X,X*,e„,e„) = -]^Hx{X,X*) = -\^{X,X*)- 

a=l 

1 ^ 1 V' 

Ri{Z, Z*) = -Y, H{Ei, E„ Z, Z*) = -Hz{Z, Z*) = -iZ, Z*), 
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and by Ri{X, Z) = 0. By assuming g = 1, we have 

(115) RiiiX) = -(^ + 0^ ; RkiZ) = -i^ + l)Z; 

(116) Rii{T) = {^+l)T, n = -{^+l){k + l + l), 

(117) Rii{X + z,x* + z*) - ^n{x + Z,X* + Z*) = 

(118) =^^ + i)^k + l-^){X + Z,X* + Z*), 

(119) Rh{X + Z,T)-^n{X + Z,T)=0, 

(120) Rh{T, T) - hz{T, T) = {^ + l){k + l + ^){T,T). 

These tensors are defined in terms of the elements of the Lie algebra. In 
order to compute them on local coordinate systems, first the metric tensor 
gij = g{dudj) , gia = g{di,da) , gap = g{da,dp) and its inverse, 9*^5*", 9"^, 
on N , need to be calculated. By the explicit form of the invariant vector 
fields we have: 

(121) Qij = 6ij + - '^{Ja{X),di){JaiX),dj), Qafj = 

a=l 

(122) = -\{Ja{X),d,), = g'^ = ^(9^, J,(X)), 

(123) g^" = 5ap + \{Ja{X), MX)) = (1 + ^x2)5„^. 

These components determine the metric tensor on SN by the formulas: 

tgij =, t^^'^gia , t'^gafd- 

7.2 Static Schrodinger and neutrino equations. 

The static model is defined by the Cartesian product, x M, of metrics, 
where M, parameterized by t, is endowed by the indefinite inner product 
{dt,dt)q = —1/q^. The several objects such as Riemann curvature can be 
computed by laws corresponding to the Cartesian products, thus they are 
non-trivial only regarding the nilpotent direction. These explicit formulas 
can easily be recover from the previous ones. 

In what follows we utilize Pauli's computation (10)-(16) regarding the 
non-relativistic approximation of the relativistic wave equation. By choosing 
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q = 1/c, the Laplacian appears in the fohowing form: 

(124) A = {Az- ^dl) + (Ax + Az + J] d^D^, • ) = 
= {Az + —dt - ^dl) + [Ax + ^x^Az + 2_, • — J-^t) ■ 

On the Z-space, operator Az — ^df^ is nothing but the wave operator (17). 
According to Yukawa's exposition, the eigenfunctions, U, of this operator 
describe the eigenstates of nuclear forces. Due to (7), the general solutions of 
this wave equation are de Broglie's wave packets (6). On the mathematical 
model, however, these wave packets are represented by twisted functions of 
the form 

(125) iPnp^,^{X,Z,t) = 

f. fc/2 

(126) =/ e'«^'^)-*)0(x,k)v.(K„)n ^1.(^)^1^ W^^ = 



(127) = / e'«^'^>-'^*V(x,k)FBp,,,(X,i^„)(iK, 



(128) where Xj^^ + 



Wave packets 'tpQpq{X, Z,t) are defined by means of the functions 

(129) Fqp,{X,Ku) = ^{Ku){%Ql){X,K^). 

They are defined also for Z-sphere bundles 5/j(x), which versions are indi- 
cated by denotations ipBp^q.SR and ipQpqSn- From this respect, ip-BpiQiR' and 
ipQpq^i correspond to the above introduced wave packets. 

For a fixed Z-vector, Z^, the regarding denotations are V'Bp,giZ^ and 
ipQpgZ^, where B is an orthonormal basis regarding Jz-yu- In this case the 
integral is taken with respect to the Dirac delta measure concentrated at 
Z^, thus these formulas can be written up without indicating this integral 

or constant ip{Ku)- By projections H^'', one defines 

(130) ^,...(X,Z,t)= /"e'«^'^>-"*V(x,k)n5^F.„(X,if„)di^, 

Jo 

where dots .... can be substituted by any of the symbols BpiC/jM', QpqM} , .. e. 
t. c, and circle, o, could symbolize any of the integral domains M', Sr, Z^. 

If also projections II^'*'' are applied to F,,, the corresponding functions 
are resp. '^^^'^^ ■ This operation makes sense only for integral domains 
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or Sji{'K) but it is not defined for the singular Dirac delta domain Zj. Anti 
de Broglie wave packets are defined by replacing —iv with +u; in the above 
formulas. The corresponding functions are denoted by ■0""** and The 
associated particles are called antiparticles. These objects can be introduced 
also by keeping —uj and replacing i by — i. 

As it is indicated, the right side of (124) is computed by adding ^^dt — 
^^dt = to the left side. Then the wave operator associated with nuclear 
forces becomes 

(131) iV = A^ + ^5,-l4, 

which is the non-relativistic wave operator established in (16). As it is 
explained in (12)-(16), the A'"-harmonic waves, defined by N {Tp){Z^t) = 0, 
relate to the relativistic wave by the formula 

(132) V'(Z, t) = e-'^^i){Z, t). 
Also remember that frequency cj is derived from 

h h \ m?c^ h h 2m ' 

by the Taylor expansion of function \/l + x. Thus the third term depends 
on h?' and, by stepping further, this exponent is increased by 2, by each 
step. For low speed particles, value u) = is a good approximation of 

the frequency, thus also E = hib \s s, good approximation for the energy of 
the particle associated with this non-relativistic wave. Note that E = E^in 
is nothing but the kinetic energy owned by the material particle. By this 
reason, the particle associated with the wave operator N can be consider as 
one of the residues of a decaying material particle which has neither mass nor 
charge and the only source of its energy is the kinetic energy of the decaying 
material particle. Such particles are the neutrinos, thus N is called neutrino 
operator accompanying the electron-positron system. 

The energy mc^ of the material particle is completely attributed to the 
other particle associated with the second operator 

(134) S = Ax + ^x^Az + ^ do^Do, • -^5*. 

incorporated into the Laplacian (124). In order to understand this particle 
represented by this operator, we introduce first the de Broglie wave packets 
^....(X, Z, t) and ^"'^**(X, Z, t) in the same way as before, but now, the iv is 
replaced by cD which can take values such as 2^k^, -^{'^f~\~'^P~\~k)pi, 2m 
4p -|- k)^ + 4/i^). By (24), the Schrodinger equation for an electron is: 
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(135) -(A(.,) - -D. . +^(.^ + y^))^ = 

On the Z-crystal model, operator <l^ is defined by the action of the Lapla- 
cian A of the nilpotent group on functions of the form tjj(X)e'^'^^^^'"^\ In 
terms of A = vrz^, this action can be described as acting only on ^ by the 
operator <^ = Ajis:+2iL'^»— /i^x^ — If the last constant term is omitted, 
the operator left is denoted by U^. Then, in the 2D case after substitution 
/i = eB/2hc, the negative of this reduced operator becomes nothing but 
the Hamilton operator standing on the left side of the above Schrodinger 
equation. If = 27rZ^, fi = k/2, m = unie, and /^(x^) is a function 
such that f^{t) is an eigenfunction of the radial Landau-Zeeman operator 
Of + 4^^, defined in (55), with eigenvalue — u) = — (4r + 4p + k)fi, then for 
5(§""|(X,Z,t)) we have: 

(136) 5(e'«^'^)+^-*)^(x2)n5')F...(X,i^.)) = 

= e^<^'^)(U, - ^|)(e^-V.(x^)<"V...(X,K„)) = 0. 

Thus on Z-crystals, operator 5* is nothing but Schrodinger's classical wave 
operator of an electron positron system. 

Note that no non-relativistic limiting was used to obtain this operator. 
It is naturally incorporated into the complete Laplacian A. Although it 
is the same as the non-relativistic wave operator obtained earlier by non- 
relativistic limiting, even the neutrino operator, N, is not the result of a 
non-relativistic limiting. The Laplacian A is the sum of these two natu- 
ral operators, meaning that it actually represents a system consisting of 
electrons positrons and electron-positron-neutrinos. The above arguments 
also suggest that this system can be regarded as the result of a sort of 
nucleus-decay. A rigorous theory describing this process is yet to be estab- 
lished. It is clear, however, that the basic mathematical tool underlying 
this physical theory must be the decomposition of the Laplacian into oper- 
ators corresponding to the constituents of a given particle system. Dealing 
with Laplacian means that one does not violates the principle of energy con- 
servation. Moreover, this tool provides also the exact operators associated 
with the particles, which is the most attractive new feature of these exact 
mathematical models. 

Actually, the elementary particles discovered in classical quantum theory 
were introduced by the very same idea. For instance, the neutrino was first 
postulated in 1930 by Wolfgang Pauli to preserve conservation of energy, 
conservation of momentum, and conservation of angular momentum in beta 
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decay the decay of a neutron into a proton, an electron and an antineutrino. 
Pauli theorized that an undetected particle was carrying away the observed 
difference between the energy, momentum, and angular momentum of the 
initial and final particles. The only difference between the two ways intro- 
ducing the neutrinos is that Pauli did not have a Riemann manifold in hand 
in the Laplacian of which he would have been able to separate the neutrino 
from the other particles resulted by the decay. 

The only term in the Laplacian containing second order derivatives re- 
garding the time variable t is the neutrino operator. This term is of first 
order in the Schrodinger operator. Because of this, waves ^ are not solu- 
tions of the neutrino operator and waves ip obtained above by the Taylor 
expansion are not solutions of the Schrodinger equation. In order to cope 
with this difficulty, non-relativistic approximation can be implemented such 
that one attributes the kinetic energy represented by Az in the neutrino op- 
erator to the Hamilton operator associated with S by considering the total 
Schrodinger operator 

(137) s = Ax + {l + \^'')Az + ^d^D^.-^du 

which is the sum of the Schrodinger operator and Az- 

In this step, the two operator is pulled together to form an operator which 
is of first order regarding the time variable. This scheme is completely anal- 
ogous to those applied by Schrodinger when, instead of the Klein-Gordon 
equation, he introduced his equation. A major difference is, however, that 
the above operator accounts also with the energy of neutrinos accompanying 
the electron-positron system, moreover, the non-relativistic approximation 
is applied to the neutrino operator and not to the electron-positron opera- 
tor. The neglected Taylor-terms in this approximation depend on h^, where 
s > 3. The wave functions regarding this pulled-together operator are de- 
fined by the eigenvalues —ij = — ((4r + Ap + k)fi + Ay?' . Then, in terms of 
FW(X,i^„) = n5'V...(X,i^„), we have: 

(138) (X,Z,t)) =S(e'(<^'-^>+2^'^*)/^(x2)F(")(X,^„)) = 
= e^<^'^')(<A - ^4)(^^^'/M(x^)F(.")(X,i?„)) = 0. 



(139) S{^^f]=S{(l e'«^'-^>+2^'^*)/ik(x2)F(")(X,K)(iK„) = 

JSr 2 

= £ e^<^'^)«., - ^|)(e^^V|k(x^)F(.")(X,i^)d/^n) = 0. 
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(140) S(§«"^i) = S( /" e'(<^'^>+2lT'^*)/ij^(x2)F(")(X,K„)(iK) = 
I e^<^'^)(<]., - ^|)(e^^*/|k(x^)F.^.")(X, i^.))k^-^rfi^„<ik = 0. 

The same formulas hold for operator S which appears as the classical 
Schrodinger operator Uij^ — behind the integral sign. Similar argu- 

ments work out also for operator (E, in which case the Schrodinger operators 
act on wave functions (x, Z, t). This is still a scalar operator which 

can be reduced to a radial operator acting on a single radial function. The 
radial operator to which the complete operator A can be reduced act on 
d-tuples of radial function, therefore integral formulas regarding these cases 
must be built up in terms of function /^II^^F^") where d-tuple (/i, . . . , f^) 
is an eigen d-tuple of the reduced radial operator. The particles defined by 
these operators are denoted by W(e = W(^(i) resp. Wa = They are 

called clean-weak and clean-high W-particles respectively, while the other 

(u) (u) 

particles = W(^(u) resp. = W^(u) are the so called dirty W- 

particles. The neutrino operator is the same in all of these cases, thus the 
associated particles are denoted Z(e. These denotations are are suggested 
by the theory of weak nuclear forces. They indicate that W-type particles 
can analogously be defined also regarding strong forces. However, the beta 
decay can be explained just by the clean weak nuclear forces. 



7.3 Expanding Schrodinger and neutrino equations. 

For the sake of simplicity the following formulas are established regarding the 
collapsing (shrinking) time-direction T under the condition q = 1- Formulas 
regarding the expanding time-direction r can be obtained by the substitution 
T = —T. Instead of t, the expanding wave functions are introduced in terms 
of . That is, the shrinking twisted wave packets are of the form 

(141) ^„„{X,Z,T)= /"e^(<^'^>"^'^'^V(x,k)nPF..(X,ir„)dii: 

Jo 

= I e'(<^'^>~'^"^V(x,k)F(")(X,i^„)dif, 



where yk^ -|- ^^^j^ = 7, and, as above, dots .... can be substituted by any of 

the symbols B^j^jM', QpqR.\ .. e. t. c, and circle, o, could symbolize any 
of the integral domains M', Sr, Z^. De Broglie's wave packets ^,,,,{X, Z, T) 
and ^""**(X, T) are introduced by the same modification, that is, the 
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ijj is replaced by u) in the latter formula, which can take values such as 
2|k2, ^{^r + 4p + fc)^*, 2|((4r + 4p + A:)^ + 4^^). 

The meson operator appears now in the form: 

(142) M = e2^Az + aT-9|y, 

The same computation implemented on the static model shows that the 
shrinking matter waves ^'....(X, T), defined in terms of are really har- 
monic, meaning M^,,.. = 0, regarding this operator. Moreover, wave packet 
^(X,Z,r) defined by 

2 

(143) ^{X, Z, T) = e-'^'^'^^iX, Z, T) 
is harmonic regarding the shrinking neutrino operator 

(144) N = e^^Az + (1 + ^e^)5T - d^^. 



According to this computation, the corresponding decomposition of the 
Laplacian into non-polarized neutrino and Schrodinger operator of a particle 
system is as follows 

(145) A = {e^^Az - 4t} + 
+{e^(Ax + ^x^Az + d^D^ * ) + + ^)^^} = 

(146) = {e^^Az + (1 + ^e^)9T - 4t} + 
+{e^(Ax + ^x^Az + Y,d^D^.) + {^ + l-l- ^e^)aT} 

(147) = (e^^Az + (1 + ^e^)aT - 9^^) + 
+e^(Ax + ^x^Az + d.D^ . -^^t) + + ^ - 1)5t. 



In terms of r = — T, these operators define the expanding non-polarized 
neutrino, Schrodinger, and tractor operators respectively. The force associ- 
ated with the third operator supplies the energy what is needed to maintain 
the expansion. Let it also be pointed out that according to these models the 
particles are not just moving away from each other but this movement is 
also accelerating. This acceleration can be computed by taking the second 
derivatives of the distance function introduced at explaining the expansion. 
This acceleration can be explained just by this new force represented by the 
third operator. 
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It is important to keep in mind that these operators are non-polarized 
The polarized operators appear behind the integral sign after these non- 
polarized operators are acting on the de Broglie waves expressed by means 
of twisted Z-Fourier transforms. 

8 Spectral isotropy. 

Spectral isotropy means that, on an arbitrary ballxball- or spherexball- 
type manifold with a fixed boundary condition, for any two unit X-vectors 
Q and Q, the Laplacian is isospectral on the invariant function spaces 3q = 

Hq ^ and Hq = satisfying the given boundary condition. Recall 

that total space Hq is everywhere dense in the straight space spanned by 
functions of the form Z){Q, X), furthermore, the boundary conditions 

can totally be controlled by {X, Z)-radial functions, therefore, this total 
function space is the same than what is defined in terms of the straight 
functions. 

Next we prove that any of the Heisenberg type groups is spectrally isotropic. 
On general 2-step nilpotent Lie groups, where the endomorphisms Jz can 
have distinct eigenvalues, this statement can be established just in a much 
weaker form not discussed in this paper. Contrary to these general cases, 
the H-type groups have the distinguishing characteristics that they repre- 
sent systems consisting identical particles and their anti-particles. Also note 
that on the expanding model this spectral isotropy explains why the radi- 
ation induced by the change of the constant magnetic field attached to the 
the spin operator is the same whichever direction it is measured from. This 
radiation isotropy, which has been measured with great accuracy, actually 
indicates that the Heisenberg type groups are enough to describe the elemen- 
tary particles and there is no need to involve more general 2-step nilpotent 
Lie groups to this new theory. 

This spectral isotropy is established by the intertwining operator lOqq : 



where heavy dot • may represent Rz{'^) or R'. They indicate the function 
spaces this operator is defined for. The very same operators are defined by 
corresponding Hankel functions obtained by the Hankel decomposition of the 



^Qpq. 



defined by 
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above functions to each other. This statement immediately fohows from the 
fact that the sums of the Hankel components restore the original functions 
in the above formulas. Also note that the {X, Z)-radial Hankel functions 
are the same regarding the two corresponded functions and they differ from 
each other just by the Hankel polynomials obtained by the projections. Thus 
the operator defined by Hankel decomposition must really be the same as 
the above operator. It follows that operator (148) preserves the Hankel 
decompositions. 

An other remarkable property of this transform is that it can be induced 
by point transformations of the form (Oggjidz), where idz is the identity 
map on the Z-space and Oqq is orthogonal transformation on the X-space, 

transforming subspace Sq, spanned by Q and all Jz„(Q), onto the similarly 
defined Sq. This part of the map is uniquely determined, whereas, between 
the complement X-spaces it can be arbitrary orthogonal transformation. 
One should keep in mind that such a point transformation pulls back a 
function just from Hg, into the function space and it is not defined for 
the whole function spaces. 

Next we prove that this operator intertwines the restrictions of the Lapla- 
cian to these invariant subspaces, term by term. Since X-radial functions 
are mapped to the same X-radial functions, furthermore, also X-spherical 
harmonics of the same degree are intertwined with each other, the state- 
ment holds for Ax- This part of the statement can be settled also by the 
formula Ax = "^dz^ idzK i written up in a coordinate system established 
by an orthonormal complex basis where Q is the first element of this basis 
{Qkui}) fo'^ ^ ^u- That is, Qq = z\ holds and the statement really follows 
by the above formula. A third proof can be derived from the fact that this 
operator is induced by the above described point transformation. 

Due to the relations 
(149) M.Fqp,(,/.) = .Fqp,((<7 - p)k</.) , Az.^Qp,(0) = .7^Qp,(-k20), 

the other parts of the Laplacians are also obviously intertwined. In these 
formulas, the second line follows from the first one by the commutativity of 
operator Dk» with the projection Hx- 

The intertwining property regarding the Dirichlet or Z-Neumann condi- 
tions on ballxball- resp. sphere x ball- type domains can also be easily es- 
tablished either by the above point transformations, or with the help of the 
Hankel transform implying that functions of the form /(|X|, |Z|) appearing 
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in the transform are intertwined with themself. Since the boundary condi- 
tions are expressed in terms of these double radial functions, this argument 
provides a second proof for the statement. A third proof can be obtained 
by the explicit formulas established for twisted functions satisfying these 
boundary conditions. 

The most interesting new feature of this spectral isotropy is that it holds 
even in cases when the space is not spatially isotropic. Let it be recalled 
that spatial-isotropy is the first assumption on the Priedmann model and the 
overwhelming evidence supporting this assumption was exactly the isotropic 
radiation measured by Penzias and Wilson, in 1965. The mathematical 
models demonstrate, however, that the spectral isotropy manifests itself 
even in much more general situations when the space is rather not spatially 
isotropic. In order to explain this situation more clearly, we describe the 
isometrics of H-type groups in more details. 

Generically speaking, these groups are not isotropic regarding the X- 
space. They satisfy this property just in very rare occasions. Starting with 
Heisenberg-type groups h'^'^\ there is a subgroup, Sp(a) x Sp(6), of isome- 
trics which acts as the identity on the Z-space and which acts transitively 
just on the X-sphere of {{^""^^'^^ . in this case the intertwining property for 
operators , also follows from the existence of isometrics transforming 

Q to Q. But the isometrics are not transitive on the X-spheres of the other 
spaces satisfying ab 0. This statement follows from the fact that the 
complete group of isometrics is (Sp(a) x Sp(6)) • S0{3), where the action of 
SO (3), described in terms of unit quaternions q by 

ag{Xi,..., Xa+b, Z) = {qXiq, . . . , qXa+tq, qZq), 

is transitive on the Z-sphere. Thus the above tool is not available to prove 
the spectral isotropy in these cases. Yet, by the above arguments, the uJqq, is 
an intertwining operator on its own right, without the help of the isometrics. 

Note that the members of a family defined by a fixed number a + b have 
the same X-space but non-isomorphic isometry groups having different di- 
mensions in general. More precisely, two members, H^'^'^^ and H^"' are 
isometrically isomorphic if and only if (a, b) = (a', b') holds up to an or- 
der. Furthermore, the sphere x sphere- type manifolds are homogeneous just 
on ~ H^'"'~^^\ while they are locally inhomogeneous, even on the 

X-spheres, on the other members of the family. Let it be emphasized again 
that this homogeneity concerns not just the homogeneity of the X-spheres 
but the whole sphere x sphere- type manifold. 
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The isometries are well known also for all H-type groups ' . The X- 

space-isotropy is obviously true also on the Heisenberg groups H^'^^ which 
can be defined as H-type groups satisfying 1 = 1. Besides this and the above 
quaternionic examples, it yet holds just on if^^'^^ ~ H'f'^\ Thus the X- 
isotropy regarding isometries is a rare property, indeed. This is why the 
spectral isotropy, yielded by any of the H-type groups, is a very surprising 
phenomenon indeed. It puts a completely new light to the radiation isotropy 
evidencing the spatial-isotropy assumed in Friedmann's model. According 
to the above theorem, the radiation isotropy seems to be evidencing all the 
new relativistic models of elementary particle systems which are built up 
in this paper by nilpotent Lie groups and their solvable extensions. These 
models are far beyond those satisfying the spatial-isotropy assumption. 

By these arguments, all the isospectrality examples established in [Szl]- 
[Sz4] for a family H^"''''^ defined by the same a + b and / can be reestablished 
almost in the same way. Note that such a family is defined on the same 
{X, Z)-space and two members defined for (a, b) resp. (a', b') are not locally 
isometric, unless (a, 6) = {a',b') upto an order. The above intertwining op- 
erator proving the spectral isotropy appears now in the following modified 
form 0,Qpq, : Sqpg, ^Qpg.' it corresponds one-pole functions hav- 

ing the same pole but which are defined by the distinct complex structures 
Jk^ resp. J'j^^ to each other. The precise correspondence is then 

(150) n^QpU^) = jS^'''^{^,K)Ilf{%Ql){X,Ku)dK, ^ 

^nr'QpM)= [ e'^^'''^ci>{^,K)uPie'^e'^)ix,K,)dK., 



which, by the same argument used for proving spectral isotropy, intertwines 
both the Laplacian and the boundary conditions on any of the ballxball- 
resp. sphere X ball- type domains. 

In order to establish the complete isospectrality, pick the same system 
B of independent vectors for both of these manifolds and, by implementing 
the obvious alterations in the previous formula, define ^Bpiqi* '■ ^BpiQi* 
Sgp.g.,. This is a well defined operator between the complete function 
spaces which follows from the theorem asserting that the twisted Z-Fourier 
transforms are bijections mapping P$b onto an everywhere dense sub- 
space of the complete straightly defined space $b- It can be defined also by 
all those maps i^Qpq» where pole Q is in the real span of the vector system 
B. Thus, by the above argument, operators ^Ispiqi* intertwine the complete 
function spaces along with the boundary conditions. 
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Interestingly enough, the isospectrahty can be estabhsh also in a new way, 
by using only the intertwining operators ^qQpq,- Indeed, suppose that the 
elements, {fp,q,i}, of the spectrum appear on a total one-pole space Sg^^ 
with multiplicity, say rupq^i. By the one-pole intertwining operators this 
spectrum is uniquely determined and each eigenvalue regarding the whole 
L^-space must be listed on this list, furthermore, the multiplicity regarding 
the whole L^-space is the multiple of these one-pole multiplicities by the 
dimension of ^Bpq- On the other hand, for Q G M''^'^", the isospectrahty 
obviously follows from Hg^g = Sg^^, thus both the elements of the spectra 
and the regarding multiplicities must be the same on these two manifolds. 

This proof clearly demonstrates that how can the spectrum "ignore" the 
isometrics. It explains also the striking examples where one of the members 
of the isospectrahty family is a homogeneous space while the others are 
locally inhomogeneous. It also demonstrates that spectral isotropy implies 
the isospectrahty. Recall that this isospectrahty is established above by an 
intertwining operator which most conspicuously exhibits the following so 
called C-symmetry principle of physics: "The laws remain the same if, in 
a system of particles, some of them are exchanged for their anti-particles." 
This intertwining operator really operates such that some of the particles 
are exchanged for their anti-particles. Thus this proof demonstrates that 
spectral isotropy implies the C-symmetry. Moreover, the isospectrahty is the 
manifestation of the physical C-symmetry. This is a physical verification of 
the isospectrahty of the examples established in [Szl]-[Sz4]. Actually these 
examples provide a rigorous mathematical proof for the C-symmetry which 
is not a theorem but one of the principles in physics. Let it be mentioned yet 
that the isospectrahty proof provided in this paper is completely new, where 
the Hankel transform appears in the very first time in these investigations. 
(All the other proofs are established by different integral transformations.) 

The isospectrahty theorem naturally extends to the solvable extensions 
endowed with positive definite invariant metrics. Just functions K) 
should be exchanged for i^(x, t, K) in the above formulas, that is, the inter- 
twining is led back to the nilpotent group. It is important to keep in mind 
that the metrics are positive definite in the spectral investigations of the 
solvable extensions. The group of isometrics acting on the sphere x sphere- 
type manifolds of SH^°''''\ where ab / 0, is (Sp(a) x Sp(6)) • S0{3), while 

it is Sp(a -|- b) • Sp(l) on SH^~^'''^\ By these formulas, the same statement 
proved for the nilpotent groups can be generalized to the solvable isospec- 
trahty family. These solvable examples provide also new striking examples 
of isospectral metrics where one of them is homogeneous while the other is 
locally inhomogeneous. By summing up, we have 
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Theorem 8.1. Operators ^Qpq and ^qQpq defined for combined spaces in- 
tertwine the Laplacians, moreover, they can he induced by point transfor- 
mations of the form {Kq^idz) resp. (Oggjidz), where idz is the identity 
map on the Z-space and the first ones are orthogonal transformations on the 
X-space, transforming subspace Sq, spanned by Q and all JzuiQ)} S'q 
resp. Sq. This part of the map is uniquely determined, whereas, between 
the complement spaces it can be arbitrary orthogonal transformation. 

By this induced map interpretation, both the Dirichlet and Z-Neumann 
conditions are also intertwined by these operators. This statement follows 
also from the fact that the very same operators are defined by correspond- 
ing the Hankel functions obtained by the Hankel decomposition of the above 
functions to each other. That is, these maps intertwine also the correspond- 
ing Hankel subspaces along with the exterior operator (E and the interior 
strong force operator S . 

So far the isospectrality is established for one-pole functions. For a global 
establishment consider a system B of k/2 independent vectors described 
earlier on the X-space. Then the global operator i^spiQi ■ T~(-^'Bpiqi{4') ~^ 
"H^gp.g. (0) can be defined also by Q -pole functions satisfying Q G Span^(B) . 
This proves that the HBptqi defines a global intertwining operator. 
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